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FOREWORD 

This report was prepared by General Mills, Inc., Minneapolis, Minnesota, 
cr. Air Force contract AP33(6l6)-7^00, under Task Nr 50855 of Project Nr 5^36, 
"Stellar Aberrascope Study Final Report". The work was administered under 
the direction of the Navigation anti Guidance Laboratory, Wright Air 
Development Division. Captain David May and Mr. Jour. Zin. provided tech¬ 
nical cognizance for the laboratory. 

The studies presented began in June i960, were concluded in December 
i960, and represent an effort of the Tracking and Control Laboratory, 
Engineering Department, General Mills, Inc. Charles Eumurian was the 
engineer responsible for directing the study. 

Although the study was a group effort, the chief contributors and their 
fields of interest were: R. Lillestrand, concept design; c. Eumurian, 

J. McGillicuday, and G. Beck, mechanization; J. Hamilton, detectors; 

R. Ginsberg, optics; R. Carlson, instrumentation; J. Carroll, error 
analysis; and C. Grosch, system analysis. 

This report concludes the work on contract AF33(6l6)-7^00, 
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ABSTRACT 


A satellite orbiting about a celestial body, equipped with a device for 
star tracking, will observe that the stars sometimes pass beh.nd or emerge from 
the limb of the body. It is shown that a measurement of the times of six inde¬ 
pendent horizon star transits will permit the determination of the six para¬ 
meters of the satellite orbit with respect to inertial space. From the orbital 
elements, the altitude and sub-satellite point can be aetermi as a function 
of time in a self-contained manner. 

Since the topographic irregularities of the planetary surface and the 
variability of the atmospheric transmission due to clouds and other sources of 
attenuation in the atmosphere might introduce large errors into the transit 
time measurement, it is proposed that star transit measurements be made while 
the light ray is still high in the atmosphei •.. The time of star transit can 
then be defined as that instant at which the refraction or the attenuation of a 
star image has built up to some preaeterained value. Each celestial body re¬ 
quires special study based on the characteristics of its atmosphere - if it has 
any. In the case of the earth an investigation has been made of the use of 
atmospheric refraction or spectral absorption at altitudes ranging from 100,000 
to 200,000 feet. As a result it appears that rms errors in the altitude of the 
constant density surface will not exceed a distance of one mile. 

In contrast wi+r. spectral absorption which is a scalar quantity, the atmos¬ 
pheric refraction is a two dimensional vector. Thus, by measuring the direction 
of the refraction in inertial space as well as using its magnitude to derive the 
transit time measurement, it is possible to derive all of the orbital elements 
with only three independent horizon star transits. If we can assume that the 
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orbit is circular tc an acceptable degree of approximation, it is possible to de¬ 
rive all o- the information necessary for a self-contained navigation system from 
orly one star transit. This can be done by measuring the two components of the 
refraction vector plus the time rate of change of refraction. 

An error analysis has shown that typically this technique will enable us to 
define our position to within about two miles for an earth satellite if each rms 
transit time error is one-tenth of a second. Error sensiti-, coefficients are 
shown for various orbital planes. Of particular note is the fact that the orbital 
period can be determined with extreme accuracy, and therefore tne semi-major 
axis of the satellite orbit can be calculated to an accuracy considerably better 
than one mile. 

By means of additional star transits or by measuring the time of passage 
into and out of the shadow caused by the sun, it is possible to define more 
than the six basic parameters. The two or three parameters of an aspherically 
shaped isorefraction surface can be derived, thereby giving an estimate of the 
shape of the central body. If the planet is a spheroid, the spatial orientation 
of the axis of symmetry can be computed. In addition, an independent determina¬ 
tion can be made of the mass of the planet. 

The refraction measuring instrument (aterrascope) consists essentially of e 
pair of opposed optical systems which are attached to a central housing contain¬ 
ing the error detection system. The identical optical systems are the Dahl- 
Kirkham type with a 100 inch focal length, 4 inch aperture and 8 inch separation 
between the secondary and the primary mirrors. Vycor brand glass will be used 
for the mirrors and 36 $ Nickel (invar) will be used for the optical system hous¬ 
ing. The wall thickness of the conical optical housing wij.l be O.Op inches, 
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this thickness will be sufficient, to meet the rigidity and vibration require¬ 
ments. The housing for the central portion of the aberrascope will be fabricated 
of lighter materials than Invar in order to reduce the over-all weight. The 
over-all length of one aberrascope is envisioned as Deing 24 inches. A system 
of two aberrascopes is pictured in the frontispiece and the o\er-all dimensions 
are 24 inches x 24 inches x 16 inches. It will probably weigh SO pounds and 
consume 80 watts of power. The measurement accuracy of the _nstrument will be 
0.1 seconds of time; it will essentially be a time measuring instrument. The 
over-all system accuracy which can be achieved is 2 miles or smaller (one 
sigma error). 

In place of rotating the entire aberrascope, the knife edge error detector 
alone will be rotated. This gives the single-axis instrument a two-axis capa¬ 
bility. The rotating knife edge can be made small and completely encapsulated 
in its owr. environment, thus avoiding the problems of rotating elements in a 
vacuum. The motor which will drive the rotor will be built right around the 
knife edge bar. The incoming stellar radiation will pass thru a window on the 
rotating error detector package and will impinge upon the rotating knife edge. 
From this knife edge, the stellar radiation will be reflected onto an ellip¬ 
soidal mirror which has alternate transparent and reflecting surfaces. Thus the 
energy will either be reflected or will pass thru the mirror. A second ellip¬ 
soidal mirror will reflect the energy which passes thru the first mirror into a 
second path. Thus the energy will ultimately arrive at one of two photomulti¬ 
plier tubes, depending on the type of surface encountered at the first mirror. 

It is difficult to abstract the 2’otating error detector: the body of this 
report must be studied in order to appreciate the salient features of this 
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device. This rotating error detector is the very heart of the aberrascope and 
its unique capability fin combination with the photomultiplier tubes) can be 
summed as follows; 

1. The incoming radiation is used 1005b of the time, 

2. One photomultiplier tube can perform the entire tracking 
function exclusive of the operation of the other tube. 

3* Chopping of the incoming radiation is perform"' 5 
rotating the image, thus simplifying the sy?• - 

h. The rotating knife-edge operates in a sealed environment. 

The entire rota+ing error detector package is translated with four servo 
drives into a null position where the stellar image impinges exactly on the 
center of the knife edge. The outputs of the aberrascope will then be the 
output of two interferometers which measure the translation of the center of 
the rotating error detector. When one star line-of-sight dips into the atmos¬ 
phere, refraction wilj build up and one end of the rotating error detector will 
shift accordingly. When this motion (refraction) has built up to a prescribed 
level (6 to 10 seconds of arc), time will be recorded aro. an occultation will 
be said to have occurred. 

The equations for determining the orbital parameters from the occultation 
times are given in the body of this report. The solution of these equations as 
well as the equations for determining the position over the earth from the 
orbital parameters will thus locate the sub-satellite point on the surface of 
the earth as well as indicate the altitude above the surface of the earth. 
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RECOMMENDATIONS 


As a result of the information contained in this report the following 
recommendations can be made at the present time: 

1) Since the potential accuracy of the aberrascope occulta- 
tion technique is about ten times greater thui the 
aberrascope aberration technique, it is recommended 
that the occultation technique be used for low altitude 
se' f-contained earth satellite navigation -yssems. 

2) Since long term thermal stability is not necessary when 
using the occultation technique, it is recommended that 
a non-rotating aberrascope design be used. This will 
result in a considerable simplification in the instru¬ 
ment design. 

3) Since the refraction measurement provides twice as much 
information at each star transit as does the spectral 
attenuation measurement and since the potential accuracy 
of the refraction technique is greater than the spectral 
attentuation measurement, it is recommended that the 
definition of star transit time be based on a purely 
refraction measuring instrument. It is also recommended 
that this instrument be designed to measure both the 
magnitude and direction of the refraction. 

U) Since the effect of instrument errors is smaller if 

we track stars to low altitudes and since at these low 
altitudes there is a concomitant increase in scattered 
background radiation from the sun, it is recommended 
that we relax our diametricity requirement in the star 
selection in order to obtain brighter stars. If we 
consider stars deviating from diametricity by as much as 
2°, we can hope to find pairs of M— 4. This reduction 
of about 3 magnitudes from the requirement for the 
aberration system will yield an increase of about 15 
times in intensity. It will also tend to relax the 
initial acquisition prob'. iras. 

5) Since the Rayleigh scattering drops off in a manner 
proportional to X ,, and since the sun's radiation 
peaks at X= p>000 A, it is recommended that we 
investigate the use of photo multipliers or photo 
conductive cells for operation in the spectral region 
from 5>000 - 10,000 angstroms. This will call for a 
search for bright low-temperature stars in contrast 
with the high-temperature stars required for the 
precise tracking when making the aberration measurement. 
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6) Since the added instrumentation for measuring the time 
of passage into the earth's shadow will be minimal, it 
is recommended that we add to the basic star tracking 
system the capability of measuring the times of ingress 
or egress of the sun, and that this measurement be made 
in the region of the COg absorption line at 4.3 microns. 

7) Since the principal value of our system lies in the fact 
that it can be made completely self-contained and since 
most viewing geometries will permit many more than six 
independent measurements, it is reccamcrde^ t^-.t vc 
immediately initiate an investigation of •••-. .-cuirements 
for a satellite-borne computer. The investigation of 
calculational procedures should take cognizance of the 
fact that the most accurate measurement will be the 
stellar transit time, the next most accurate the solar 
transit time, the next most accurate the direction of 
the refraction vector in inertial space and the least 
most accurate will be the measurement of refraction rale. 

8) The measurement of the magnitude of the refraction to an 
accuracy of 0.25 to 0.5 seconds of arc will be sufficient 
to keep the instrument error equal to or less than the 
meteorological errors. However, the measurement of the 
direction of refraction should ideally be good to about 
1/2000 radian in order to yield an error of 2 miles on 
the earth's surface. This implies an rms accuracy of 
0.02 seconds of arc in a 40 second of arc refraction 
measurement. Because of the desirability of using the 
direction as well as the magnitude refraction information, 
it is therefore recommended that we attempt to achieve 

an accuracy 01 0.05 to 0.1 seconds of arc when measuring 
the atmospheric refraction with the non-rotating aberrascope. 
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INTRODUCTION 

This report is tne second and final report being submitted to the Wright 
h-i Development Division in accordance with the requirements of Air Force 
contract rumber AF 33(6d6)-7^00 which is entitled: "Theoretical Investiga¬ 
tion of New or Unexplored Techniques for Recoverable Orbital Vehicles". 

The :revi.as quarterly report presented an analysis of the guidance problem 
in the light of an aberration sensing instrument. The of the 

various aspects of an aberration measuring instrument were given. At the 
same time, a preliminary analy:is was given to the consideration of an 
instrument which measured occultaticn time for a given star. From the 
cccultation time measurements, the orbital parameters and finally the 
position over the earth can be established. 

While It was stated in the last quarterly report that the instrument 
design under consideration could be used for either the aberration or the 
occultation techniques, further study has shown the feasibility of chang¬ 
ing some of + he original instrument design characteristics in order to take 
full advantage of the occultation technique. This technique does have dis¬ 
tinct advantages and its use will simplify the instrument design but will 
limit its use to the occultation method. The body of this report describes 
in detail the phenomenological basis on which the occultation measurements 
can be made. In addition, the salient features of the occultation instru¬ 
ment will be pointed out. The instrument is now called the "non-rotating 
aberrascope". In general two of these units, mounted 90° to each other, 
will be used for determination of the orbital parameters. 

Again, as predicted in the last quarterly report, this report 
emphasizes the actual design of an instrument. Drawings of various 

portions of the instrument as well as detail specifications are 

"Manuscript released by the authors, December i960, for publications as a 
WARD Technical Report." 
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included in this report. Some Quotations on the cost of hustings and 
other out-of-house purchased parts have been received, i'nus, we are 
well prepared to jove ahead into the prototype development as the 
next logical step in this satellite guidance program. 
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II. GUIDANCE SYSTEM BASED ON MEASUREMENT OF STELLAR OCCULTATION TIME 

A non-rotating aberrascope can be used to measure the times of transit 
of a star behind the earth, and it can provide information regarding 
the direction of the refraction in inertial space. Nith these measure¬ 
ments it is possible to determine all of the orbital elements by observing 
three star transits. If we add a number of simple infrared sensitive 
phoi-oconductive cells to measure the time at which r itellite 
passes into and out of the earth's shadow, only two star transits are 
needed to determine the orbital elements. 

By measuring the refraction rate or the intensity loss due to 
differential refraction, additional information can be obtained 
regarding the altitude of the satellite. These two measurements lead 
to greater errors than the stellar or solar transit time measurements, 
so they are expected to be used only in making first approximations 
of the orbital parameters. 

A. Transit Time Measurement Only 

An observer stationed on a satellite can determine his orbital 
parameters by observing the times at which at least six known stars are 
occulted by the planet about which he is rotating. Having determined 
the parameters of his orbit he can predict his position relative to the 
planet as a function of time, and he therefore has a completely self- 
contained navigational system. The calculations and discussion given 
below reveal some of the factors involved in designing such a system. 

The first problem encountered with a navigation system based on 
the measurements of stellar occultation times is that of the irregular 
shape of the earth and the presence of clouds at various altitudes. 

These effects are illustrated in Figure 1. Since the highest mountains 
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rise co more than five miles, a limitation to ac u>~acy is immediately 
imposed. Clouds present an even greater source of error in that they 
are sometimes observed at altitudes in excess of ten miles. These 
two sources of error could be circumvented by using for occultations 
some phenomenon which lies at altitudes in excess of those where 
clouds are observed. 

Two possibilities are shown in Figure ±. One >„■ these is based 
on the measurement of atmospheric refraction and the other on the 
measurement of the atmospheric absorption of starlight in various spectral 
regions. These techniques will define an "occultation surface" surround¬ 
ing the earth. This surface will be defined as the point at which the 
atmospheric refraction or the spectral absorption builds up to some 
predetermined value. 

To obtain the ultimate accuracy achievable with the navigation 
system, an occultation surface must be selected possessing an altitude 
characteristic which can be accurately predicted. For example, if an 
occultation surface is used based on the absorption caused by ozone 
in the .20 to .29 micron region of the spectrum, an attenuation of the 
signal from the star of 50$ might, on the average, be observed at an 
altitude of 115,000 feet (35 km.). But it is well known tbat the 
ozone composition of atmosphere undergoes extreme latitudinal and 
seasonal variation, and in order to schieve the best possible navigational 
accuracy it would be necessary to correct for these var.ations. If the 
altitude characteristic of the "occultation surface" were defined on 
the basis of atmospheric refract;on, corrections due to seasonal and 
latitudinal variations in the temperature and density of the atmosphere 
would also have to be made; however, the isorefraction surfaces are not 
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expe:~ea ~c fluctuate as widely as the ozone isoabsorption surfaces. 

It the satellite were rotating around a planet with an unknown mass, 
and if the altitude of the occuitation surface is not known, an additional 
^.-.'altaticr. measurement would be required for each of the additional para- 
me ’c.n tc c. determined. For example, if one assumed that the occuitation 
surface hai *ne shape cf an ellipsoid of revolution, two additional occulta- 
tious would be needed to define the parameters of this surface. If the 
pj.ar.et itself nad an aspneric shape, this could also be determined with 
additional stars, however, a mere accurate determination could probably 
to strained by noting the secular variations in the longitude of the node 
and in + he argument of the perigee over a large number of orbital 
rotations 

Since the semi-major axes of the orbit is a function of the oreital 
period only, we have 


a = Semi Major Axis 
T = Period of Orbit 
fe = Radii.t. of Earth ^ 

g = Gravity at Earth's Surface 

and therefcie by measuring the time from one occuitation of a given star 
tc the next occuitation of this star, we can computer the semi-major 
axes m a particularly simple manner. In addition, the systematic 
errors in the measurement of occuitation time cancel out in determining 
the orbital period, and thus this parameter can be determined with great 
accuracy. 

Frcm Equation ( i , we can determine the accuracy with which the occulta- 
tion time must be measured. If g(a) is the rms error in the calculation 
cf the semi-major axis an> o' (t) the rms error in the measured orbital 
period, we can write 

3"(T) _ 3 g (a) 

T 2 
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For low altitude earth satellite orbits, a ^4,000 miles and T ^ 100 
minutes. Since the errors caused by the uncertainty in the altitude 
of the occultation surface are of the order of l/2 mile, the accuracy 
with which Qr (T) should be measured in order to contribute a 
comparable error is 0~ (T) 1 second. 

The occultation technique finds its greatest value for low 
altitude satellite orbits. This becomes evident fr F_g*.re (2) 
where we have shown the relation between the orbital altitude and the 
percentage of the sky occulted by the figure of the earth during one 
rotation of the satellite. For a circular orbit of 400 n. miles 
altitude, stars lying in 90 per cent of the sky will be occulted 
during each rotation, while for orbits of 40,000 miles altitude,only 
about 8 per cent of the sky is swept out during each rotation. This 
altitude problem is further compounded by the fact that orbits which 
are inclined relative to the earth's equator will undergo a nodal 
regression due to the oblateness of the earth. This results in 
certain periods in which those stars whose position., are close to 
the ay.is of the orbital plane will undergo no occultation. 

B. Maximum Information Derivation 

Thus far we have considered the use of the occultation technique 
to form a self-contained navigation system in which only occultation 
time measurements are used to define the orbital parameters. The 
refraction measurement provides no more information than the spectral 
absorption as a means of defining the occultation event. 

Refraction may, however, be regarded as a two dimensional vector 
having direction as well as magnitude - in contrast with spectral 
absorption which is purely a scalar quantity. Thus, the following 
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information can be derived from each occultation measurement if the 
atmospheric refraction is measured. 

1) Occultation time, seconds, T 

2) Direction of refraction, radians,^ 

3) Time function of magnitude of refraction, (\ (t) 

These are listed in order of accuracy from the standpoint of the self- 
contained determination of the orbital elements. If we are dealing with 
a circular orbit , it is possible to determine all Cl u:.e elements of the 
orbit from a single star transit measurement in the following manner. 

1. Time Measurement 

As a result of the measurement of the star transit time, we 
know that we are somewhere on the surface of an elliptical cylinder 
at the time of the star transit. This is shown below. 



The direction of the axis c this cylinder is given by the position 
of the star we are tracking, its shape, by the shape of the isorefraction 
layer and, ultimately, by the shape of the earth. 

The errors in the time measurement will not introduce an error 
of more than one or two miles. 

2. Direction of Refraction 

Having established that we are somewhere on a certain cylinder, 
the measurement of the direction of the refraction enables us to deter¬ 
mine that we are somewhere along a line on this cylinder - this line being 
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a straight one lying parallel to the axis of the elliptic cylinder. 



If we can measure the refraction to O'.'l when the magnitude is 50'.'0, 
then the error in direction will be of the order of 1 /, '0U reins. At the 
radius of the earth this will correspond to about 3 C XJ/500 = 7 n. miles 
error in position. 

3. Time Function of Magnitude of Refraction 

By measuring the rate of change _f refraction it is possible 
to estimate the sem.. ^ajor axis of the orbit (assumed in this case to be 
circular). 

Z W 2 ) 1/2 



S 
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as follows. 

The atmospheric refraction may be written 
R =K 0 exp[-K,(/> t - /?>] 

R =-K 0 RR t = .K 0 R AT 0 

••-v o =-R/k q R 

Since we can measure R and R , we can compute / V Q . Knowing 
3T we can compute /& . 
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At R - 50 sec., the value of /? is probably of the order of 
50 sec/sec. Probably £R could be held to 0.1 sec., but<$~/? will 
probably be of the order of 1 sic/sec. 

This means that - _ — ^ 2 x 10 . From Figure (4) it can be shown 


V 0 

that for a 300 mile orbit ^ -AT 

Z h 

Since AT Q = 9,000 ft/sec 
<£/ir 0 = 180 ft/sec 


2,000 ft/'sec = 17 ft/sec/mile 
150 miles 


and therefore £ h = 1 mile (l 80 ft/sec) 

17 ft/sec 

h i 10 miles! 


C. Combinations or Measurement Techniques 

The basic navigation technique under consideration i3 that of the 
measurement of star transit times, but the foregoing discussion leads 
to the consideration of several other measurements which, when combined 
with the basic technique, yield either 

1) additional parameters, such as the mass of the planet or 
the shape of the .transit surface, 

2) a complete determination of the six orbital parameters in 
less than one complete rotation of the satellite, or 

3 ) a more accurate determination of the orbital elements by 
selecting the calculation technique yielding the smallest 
values of the error sensitivity coefficients. 

Figure ( 3 ) contains a summary of the various combinations of 

measurements. In this figure one series contains the cases in which 

only the times of ingress are measured; another series contains the 

cases in which the ingress times plus the times of egress occur in the 

earth's shadow (reacquisition is a problem if the egress occurs in an 

atmosphere illuminated by the sun); the last case contains the general 

solution in which all ingresses and egresses are included. In each of 
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these cases there may be eltner tvo or four stars occulted per rotation. 
This occurs because the stars are picked in diametrically opposite pairs. 

In o' der of accuracy for instrumentation now under consideration, 
the following additional measurements can be made: 

1) time of passage into earth's shadow, 

2) direction of i" fraction in inertial space, and 

3) refraction rate as a function of time- 

Various combinations of these techniques are shown in Figure (3)• 

For example, if we measure the time of ingress only, if four stars are 
occulted per rotation and if we measure the direction of the refraction 
vector at each ingress, then Figure (3) shows that nine measurements are 
available per satellite rotation. 

The use of the refraction rate information is n~t shown in Figure (3), 
since this information is not very accurate and will be used primarily for 
the first estimate of the orbital parameters rate- r than for the final 
calculation. In deriving Figure (3) it was assumed that at the completion 
of 360°rotaticn one additional measurement would be used to determine 
the semi-major axis. It is not known whether all of the combinations 
shown in Figure (3) are independent and therefore additional analysis 
is needed regarding this question. 

The measurement of star transit times alone leads to certain 
conditions in which che orbital parameters either cannot be determined 
or are subject to large errors. Thus, pernaps the most important point 
to be made as a result of Figure (3) is that the addition of refraction 
direction and solar transit time measurements results ir. an extremely 
flexible system with many alternative methods of calculating the 
orbital parameters. The price of this flexibility is very small, in 
that the additional measurements require little additional instrumentation. 
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III. 


PHENOMENOLOGICAL CONSIDERATIONS 


A. Methods of Defining Transit Time 

As starlight moves into the earth's atmosphere or emerges from 
behind the earth, the dynamic response of the detection instrument 
must not limit the measurement of time. The faster the rate of 
ingress or erress, the more difficult the instrument design problem. 

Consider the geometry shown below, where the ing symbols 

radius of planet 

radius to occulting layer 

radius to satellite 

angle from point of tangency of ray 
on occulting surface to satellite 

gravitational acceleration at^ 

satellite altitude = 

velocity of satellite 

component of satellite velocity 
perpendicular to occulting surface 
at point of tangency of ray 

For the present example, assume that we are dealing with a 
spherically shaped planet and a circular satellite orbit. 



Ilf 


are used: 


/O _ 
t o - 

a - 

e = 


rr 






^ =K.sin 9 


but sin 9 = () 


1/2 




1/2 4 


Therefore 

_ 1/2 2 j2 1/2 

^ O ■/_£ : .8_l/ ° s ~_/ ) - 

A m 

To find the maximum value of AS' ve let <)/«■ o =0, 

o 

And therefore - ~\J3^ ■ For the earth = 3^50 ■} 30 

n. miles, and therefore the velocity /V will be a maximum for 
satellite orbits of /° g = 6030 n. miles, or satellite altitude = 2580 
n. miles, above and below this altitude /tf will be smaller, and there¬ 
fore the instrument dynamic response problem will be easier. At 
h = 2580 n. miles, the value of (/U~ 0 ) max. = ]6,000 feet/sec. The 
value of AT at various satellite altitudes is shown in figure (h). 

The nature of the function is such that /(/• Q varies surprisingly 
little over the usual range of altitudes considered for earth 
satellites. 

I. Atmospheric Refraction 

One of the modes of operation of the system previously 
described is that of defining the star transit time on the basis of 
some predetermined magnitude of atmospheric refraction. As will be 
eviuent from the charts shown later, the change in the refraction 
caused by the atmosphere is so fast compared with the change in the 
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Figure 






ve]■'city aberration -chat the latter may be neglected. 

Figure ( 5 ) shows the approximaxe refraction for rays grazing 
through the atmosphere at various altitudes. At an altitude of 
200,000 feet the refraction was assumed to be 1 second of arc, while 
at 100,000 it was assumed to be 60 seconds of arc. If we assume an 
exponential atmosphere, the relation between refraction and altitude 
can therefore be written 

ft = 3-64 x ic3 exp (-4.10 x 10"5 h) 
where ft is in seconds of arc and in feet. 

In selecting the magnitude of the refraction, ft Q for the 
transit time measurement one must take account of the following 
factors. 

1) How predictable is the isorefraction surface in question? 

2) Are extensive seasonal and altitudinal corrections necessary? 

3 ) Are these corrections known? 

4) How accurately can the star tracking instrument measure 
the refraction? 

5) What tracking problems arise because of the scattered 
radiation from the sun? 

6) Is the atmospheric attenuation of the starlight sufficiently 
small at the level in question to permit ft to be measured? 

Requirements ( 5 ) and ( 6 ) can be best satisfied by selecting an 
occultation altitude which is extremely high. Requirements (l), ( 2 ), and 
( 3 ) tend xo disqualify both the very high and the very low altitudes. 
Requirement (4) disqualifies the very high altitudes. Thus, if we 
are to make the optimum selection of the altitude of the occultation 
surface we must pick one lying at an intermediate altitude. On the 
basis of the factors mentioned above, it would appear that this surface 
should lie between 100,000 and 200,000 feet above sea level. 
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To illustrate how the accuracy of the star tracking instrument 
effects the altitude measurement, we can derive the refraction- 
altitude sensitivity coefficient from the above equation. We then 
obtain 

ah = - 6-70 exp (4.1 x 10-5 h) f e et 
~w sec. of arc 

This sensitivity coefficient gives the number of feet of errc v ' in the 

occultation surface resulting from a given error j.n Lut measurement 

of the atmospheric reiraction. 

Figure (6) is a graph of this equation. It shows that an error 
of 1 second of arc causes an altitude errur of only 400 feet at an 
altitude of 100,000 feet. As the altitude increases, the effect of 
an angular error in the measurement of the refraction increases very 
fast - particularly above 150,000 feet. Thus unless an extremely 
accurate instrument for angular measurement is used, it would be 
desirable to use an occultation surface lying below 150,000 feet. 

On the other hand, if an occultation surface lying below 150,000 feet 

is selected, the response of the instrument must be very fast. This 

can be seen from Figure ( 7 ) where the magnitude of the atmospheric 

refraction as a function of time is plotted on the left hand ordinate. 

On the right hand ordinate the refraction rate is plotted as a function 

of time. This is derived from -the equation 

d = - .149 / dh ) exp (-4.10 x 10^ h) seconds arc 
dt ^ dt ' second time 

in the case of Figure ( 7 ) an orbit of 300 n. miles altitude and an 

occultation surface of 150,000 feet were assumed. 

Had an occultation surface lying at 130,000 feet been selected, 
it would be necessary to make the refraction measurement while it was 
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changing at a rate of about 6.3 seconds of arc per second ,-f time. This 
raises a question of the possibility of building an instrument which 
measures the rate of change of atmospheric refraction rather than the 
absolute value of atmospheric refraction. 

The remainder of the problems associated with the evaluation of the 
six factors affecting the selection of the altitude of the occultation 
are meteorological in nature and are diseased in :: -uesequent 
paragraph. 

2. Atmospheric Attenuation 

It is possible to define star transit on the basis of the 
spectral intensity rather than the atmospheric refraction. Some of 
the problems associated with this measurement are considered in 
this paragraph. 

First, it should be pointed out that the spectral intensity is a 
scalar quantity while atmospheric refraction is a vector quantity. Thus, 
if we measure spectral intensity at each transit we obtain only half 
as much information as we do by measuring atmospheric refraction. 

The received signal at the telescope will vary in intensity at a 
given altitude of transmission through the atmosphere for a variety of 
reasons. Some of these are listed below: 

1) absorption of stellar radiation by atmosphere, 

2 ) scattering of stellar radiation by atmosphere, 

3 ) differential refraction of stellar radiation by atmosphere, 

>0 scattering of solar radiation by atmosphere to create 
background noise, 

5 ) radiation from the atmosphere itself, 







6) large scale variation of absorbing and scattering com¬ 
ponents such as Oj, HgO, and small particles, 

7) small scale variation of meteorological conditions which, 
to some extent, can be avoided by seasonal and latitudinal 
corrections. 

In general, all of these phenomena can contribute errors to the 
transit time measurement. The following paragraph discusses the 
selection of a spectral region and occultation altitude which will 
minimize these errors. 

a. Absorption 

We would like to pick a spectral region in which some¬ 
thing like 50 per cent absorption* occurs at an altitude in excess 
of 60,000 feet in order to make the definition of star transit time 
independent of .local cloud cover and meteorological variations. This 
requirement immediately restricts us to the ultraviolet region of the 
spectrum. In the visual portion of the spectrum the absorption is 
too small; in the near infrared, where 1^0 or COg absorption bands 
might be used, the intensity of starlight is not adequate. 

In the ultraviolet region of the spectrum the absorption resulting 
from Og and 0^ is shown graphically in Figures (8) and 9)- 

(l) Ozone - A gaussian distribution of ozone was 
assumed having its peak at 25 km and a JL width of 10 km. Ther, 
using the absorption coefficients for A = 2,000-3^000 A, 1 the per 
cent transmission was plotted against h with A a parameter. For 
this wave length range, we find that .25 JJ has the cutoff highest in 

* We could consider the altitude at which 50 per cent scattering takes 
place (or a combination of scattering and absorption), but if the 
speci-ra! region in question gives a large scattering of starlight, it 
also yields a large background intensity due to scattered sunlight. 

Tnis is not desirable. 
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tne atmosphere. This is snown in figure (8). 

(2) Op - Again using the air masses calculated, in the sec¬ 
tion on scattering, and assuming that 0 Q constituted 22 per cent of this, 
and using absorption coefficients for 0^, the percentage transmission is 

O 

plotted against h, against using X as a parameter. Here X = l lf -5 0 A is 
the wave length yielding the maximum cutoff altitude. This is shown in 
Figure (9)- 

From Figures (8) and (9) the following summary is made of the alti¬ 
tudes at which a 50 per cent atmospheric absorption for radiation in the 
various regions of the spectrum occur. 


O 


SPECTRAL REGION, A 

ALTITUDE, KM 

PRIMARY ABSORBING 
CONSTITUENT 

1330 - 1650 

72 

°2 

1360 - 1540 

79 

°2 

1450 

82 

°2 

2200 

46.5 

°3 

2500 

50 

°3 

2800 

48 



In general, the altitude for 5° per cent absorption increased as the 

O 

near ultraviolet at 3 j 000 A is approached. Were it not for the extreme 
fluctuation in the atmospheric ozone content, the region from 2,200-2,800 
angs'.roms could be very attractive. While the 0^ content does not fluct¬ 
uate to an objectionable extent, its primary absorption lies in a region 
of the spectrum which is largely unexplored. A star would require a tem- 

O 

peradure of 20,000° K to peak at 1,500 A, and it woula be desirable to use 
0 or B spectral classes if this region of the spectrum were used. For the 
stars down to 8.5 visual magnitude, approximately 10 per cent are 0 or B 
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stars, and thus the selection of stars is considerably restricted by 
operation in the 0^ regions of the spectrum.^ 
b. Sc attering 

Atmospheric scattering has two effects. It creates 
background radiation and it attenuates the signal being detected. The 
following paragraph discusses the latter effect. 

Rayleigh's scattering law indicates xnat the .si.sity of scattered 

\ I* 

radiation is proportional to A “ r - Thus the region in which the 
highest absorption levels are obtained at the short wave lengths is 
also the region in wh.ch the greatest scattering attenuation is 
obtained. However, as is evident from Figure (].C,, the signal 
loss due to atmospheric scattering is not appreciable above 30 tan. Since 
absorption is 99 per cent for Og by 60 km and 0^ by bo km, it is 
evident that the scattering attenuation may be neglected. 

In deriving Figure (10) observations were made which deserve the 


following comments. The number of air masses encountered in a 
horizontal line through the atmosphere at a height h was calculated. 


Then using P. Moon's^ results 


sp, , (bO 

T- (f a A> =10 



This gives the transmission coefficient, I^ , for scattering as a 


function of the scattering coefficient, 7^/1 > normal pressure 


and temperature for one air mass, m is the number of air masses 


calculated above, p is pressure in mm at height h, ana is taken 
as a function of wave length. 7 

Since p end m are functions of h, and k is a function of /? , 

we plotted in per cent against h using A as a parameter. 

How p is really too high since it is calculated at the point of 
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greatest density along the path. But for smaller p, the graphs 
merely get moved to lower heights. 


c. Differential Refraction 

Because of the fact that refraction increases as a ray 
dips into the atmosphere, the intensity of the ray is reduced. The 
atmosphere acts like a negative lens and disperses the radiation. 

Under certain conditions the loss in signal strengc . uae zo differential 
refraction is greater than that resulting from absorption and scattering. 
In general, the effect of differential refraction becomes greater as the 
altitude of the satellite increases. 

In order to show how this occurs, consider the sketch shown below. 



Elemental area before 

refraction, dA 
' o 



Elemental area after 
refraction, dA^ 



The change in intensity is equal to the ratio of the elemental areas 


. Thus I = dA 0 


dA n 


dA, 


ZjLi£jL£_-= 

y dy d 9 ' y ay ' 


(1 - D R ) ( 1 - D d'/? ) 

~T dp~ 

But f^= K q exp (-K-^ ~ {O 7 ) and therefore 

_d_ P- = - K. @ and I = _1_ 

d J~ Q 1 (1 - D /P~H 1 + /? ) 

f c 
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Since we are not interested in satellite altitudes greater than ICr miles 
or in atmospheric refractions greater than 100 seconds of arc, the maximum 


value of D_ is approximately 


f S x 

104,000 x 100 x 4.85 x 10"° 


- 1.46 x 10' 2 


We may therefore use the following approximation 




The values of f? and d R can be obtained from eolations given earlier. 
d /o 

’/hen expressed in units of radians and nautical miles we may write 

^ = 1.765 x lCf 2 exp [ - .249 {j? Q - j' e )J radians 
d i? _ -.249]? radians 

Tfa “ n - mile 


V en these values are graphed. Figure (ll) is obtained. This figure 
shows that if the ray penetrates to 100,000 feet (corresponding to about 
60 seconds of arc refraction), the effect of differential refraction can 
be significant - even for low altitude orbits. Using the ozone absorption 
with 50 per cent attenuation at about 160,000 feet, a satellite altitude of 
1000 n. miles will produce an effect of differential refraction not exceeding 
2 per cent. 

A slightly different form of Ihe in ensity function is shown plotted in 
Figure (12). A satellite altitude of 30C milts was assumed and the rela¬ 
tion between height and per cent transmission was derived, 
d. Scattered Radiation from the Sun 

When the viewing geometry is such that ingress or egress takes 
place in an atmosphere illuminated by the sun, the background radiation must 
be considered. Scattered radiation from the sun presents the greatest pro¬ 
blem in the visual region of the spectrum where the maximum of the solar radia¬ 
tion lies. For wave lengths below the solar maxima at,A = 5>0CO A, the inten¬ 
sity per unit wave length t3y be approximated by the Wien radiation formula 
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which states I 




Kn 


A ' exp ( ch_) 

n Jk T 1 


The Rayleigh scattering law states that 




A 

Therefore, since net scattering is proportional to the amount of 
solar radiation to be scattered times the scattering coefficient, we 
obtain as a measure of the scatter intensity 


= I 


■<A W = 


K- 


1 


/\' exp / ch 


G-w-) 


T'-. find the wave length at which the maximum cf I is reached we let 




= 0, and obtain the equation ch = 9 

AMt 


Since ch = 1.1*39 x 10^ when A is in microns and T in ° K, we have 
k 

for the sun, letting T = 6,000° K, 

A Q = -2665 microns 

Thus, a wave length of 2,600 A presents the least favorable 
situation regarding scattered radiation from the sun. Of course, 
this conclusion applies per unit wave length only, and it may be 
expected that the effect of scattered solar radiation will be 
greatest, i.e. least favorable, in the ozone region of the spectrum. 

At a = 1,500 A in the Og region the solar effect per unit wave 
length interval will be reduced to 


I ( ?[ = 1,500) / -266 \ 9 ( l 9 ' 00 ) 176 

I ( A ~ 2 , 660 / “ t” 150 ) IplSTCOv- 11,000 


1.6 X 10 


-2 


Thus, the Op region has a considerable advantage over the 0^ region. 
Figure (13) shows the effect of scattered solar radiation for a 

O 

spectral region from roughly 2,500 to It,000 A. The FOV of the 
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telescope was assumed to have a diameter of only 5 seconds of arc - 
a value which is probably the best that can be achieved. Since 
a star of B spectral class and 3rd magnitude will yield approxi¬ 
mately 10^ photoelectrons per second, it is apparent that scattered 
radiation will become troublesome at altitudes below 163,000 feet (as 
many photoelectrons as in the signal). Since the 50 per cent 
absorption point for ozone occurs at about 160,000 feet, it is 
apparent that solar scattering will be a serious problem for 
operation in the ozone region of the spectrum. 

e. Emission from the Earth's Atmosphere 

Emission from the earth's atmosphere in the wave 
length regions of interest is so small that it can be neglected. 

At very high altitudes in the thermosphere, high kinetic temp¬ 
eratures are observed ; however, the atmospheric density is so small 
that the radiation is insignificant. At the lowest altitude, 
the ray will pass (for ozone) at the level of the mesopause where 
the temperature reaches a maximum of about 300° k. At this temper¬ 
ature and for /[ 2 microns, only 10“5 of the total self-emission 

is included. 

f. Large-Scale Variation of Absorbing Components 

Of the gases of the atmosphere which undergo large- 
scale variations, 0^ and HgO are the most important. Generally, the 
variation of 0^ is sufficiently great to result in rms errors in the 
altitude of the isoabsorption level of at least 3 or *4 miles. The 
major HgO absorption bands lie in the infrared and have not been 
considered as a source of stellar attenuation. 

Additional data describing the atmospheric ozone fluctuation is 
contained in a later paragraph. 
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g. Small-Scale Variations of Absorbing Components 


The other gases of the atmosphere (other than 0^ and 

H 2 O) are not subject to such extreme variations and thus their 

use as absorption media is considerably more attractive. The major 

factor affecting the absorptivity of these gases are the diurnal, 

seasonal, and latitudinal variations of density at any given level. 

B. Msterological Data 

1. Atmospheric Density and Temperature 

The classical s'.eory of refraction for a ray passing 

through a planetary atmosphere provides the expression 

t h 

R = ( n c - 1 ) I" 27 r £0 1 ' 


rim 

where H = _ 


where R = refraction in radians 

-1 = modulus of refractivity at the altitude 
j> Q - distance of closest approach of ray to center of earth 
R = gas constant 
T = absolute temperature 

M = mean molecular weight of atmospheric gases 
g = gravitational acceleration at j> Q 
Since (tt Q -l) is proportional to the air density, D, and since 
the mean molecular weight M may be assumed to be constant at the 
altitudes in question, the atmospheric refraction is 

TO K ° D 

^ 172 

T 

Meteorological data indicates that at any given altitude the 
statistical fluctuation in temperature is smaller than the 
statistical fluctuation in density. If we assume that the 






statistical fluctuations in D and T are not correlated we may write 





1/2 


Because CT (T)/T is smaller than <T (D)/D, and because T appears 
under the square root thus generating the factor l/4, we may 
estimate the rms fluctuation of the refraction from the rms fluctu¬ 
ation of the density as follows: 

^] 

The problem of determining the rms error in the occul-ca-cion altitude 
now can be expressed as a problem in determining the statistical 
fluctuation in the altitude of the isopycnic surfaces (D = const.). 

In making estimates of the atmospheric refraction for a grazing 
ray, the first approximation can be obtained by multiplying the 
sea level refraction at 90 ° zenith angle by twice the density ratic. 
/?(h) » 2 fl(0) ^_Do 'j 


<T ( !< ) = R 


[■ 


This has been done in deriving the curves shown earlier. 

The next refinement in the refraction calculation is to correct 
for the scale height, H, at the various altitudes. Since H = RT/Fig, 
we must consider the variation of T, g and M with altitude. Below 
100 km we may assume that M = 28.97* 50 that the principal sources 
of variation lie in g and T. 

Newall^ gives the values shown on the table on the next page. 
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IIGHT (km) 

TEMPERATURE (°K) 

g(cm/sec 2 ) 

SCALE HEIGHT (km) 

1.216 

291.0 


8.53 

c: 

J 

276.8 


7.83 

10 

230.8 


6.78 


209.I 


6.l6 

20 

212.8 


6.28 

25 

223.0 


6.59 

30 

231-7 

97 c.4 

6.85 

35 

244-5 


7-24 

ho 

262.5 

967.3 

7-79 

45 

271.3 


8.06 

50 

270.8 

964.3 

8.06 

55 

265.8 


7.93 

6o 

252.8 

961.3 

7-55 

65 

235.0 


7.03 

70 

218.0 

958.4 

6-53 

75 

209.1 


6.27 

80 

205.O 

955-4 

6.16 


Since the refraction R> is proportional to the inverse square root of H, 
an inspection of the above data shows that the maximum error in 


using the assumption of constant H is less than 20 per cent. This 
would correspond to an altitude error of about 5>000 feet at the 
altitudes where we will be working. 

Since the refraction of a grazing ray is related to the vertical 
density gradient, the fact that the atmosphere is not isothermal creates 
an additional complication because the existence of a /' c) T ) 

TK 

term creates a source of density gradient other than the exponential 


38 




atmosphere. In the final case it will be necessary to integrate 
numerically over various atmospheric paths. 

Figure (]b) shows the data contained in Figure (6) expressed in 
terms of the error in feet resulting from various errors in the 
measurement of atmospheric refraction. Also shown in Figure (l4) 
are the approximate errors in the altitude of the isorefraction level 
due to meteorclogical fluctuations. The two extremes represent 
the case in which the best available latitudinal ce'^unal corrections 
art used and the case in which no seasonal or latitudinal corrections 
are used. These limits have been estimated and are described later in 
this section. 

It is evident from Figure (l4) that an instrument error in the 
measurement of atmospheric refraction of 1 to 2 seconds of arc at an 
altitude of 100 to 120,000 feet will lead to an error which is 
comparable with the meteorological ej rors. However, there are several 
reasons why operation at these altitudes is not desirable. They are 

1) As shown in Figure ( 13 ), background radiation from the 
sun is very great. Probably only measurements in the 
sun's shadow can be made at these altitudes. 

2) At 100,000 feet the rate of change of refraction is 
about 22 sec of arc/sec time, and the dynamic response 
of the instrument must be very rapid. This requires 

a wide band width; however, a wide baDd width is 
antithetic to precise tracking in the presence of 
large amounts of background noise 

3) Since we wish to measure both the time of egress as well 
as ingress, we must be able to acquire and stably track 
an emerging stsr at 100,000 feet. This will be extremely 
difficult if we have a small FOV (necessary for tracking 
in the presence of large amounts of background radiation). 

The net result of these factors is that the transit time should 
be defined for altitudes in the neighborhood of 150,000 feet. If we 
wish oO keep the instrument errors small compared with the meteorological 
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errors, and if we make nominal seasonal and latitudinal corrections, 
it is apparent from Figure (l4) that an instrument with an accuracy 
of 0.25 to O.p seconds of arc can be justified. 

a. Height Variation of 30 km Surface 

This paragraph discusses the problems of fir-ding the 
seasonal and latitudinal height variation of a density surface which 
varies about the 30 km level. The data selected for use in this 
-tudy is published in the AB4A Climatological Rir.gboo-' ‘ This 
report gives monthly table., of mean and standard deviations of 
density for various altitudes over ten selected stations. 

The monthly values listed at the 28 through 31 km levels were 
combined into seasonal values. This was done for only nine stations, 
since one of the ten did not have sufficient values at these levels. 
The mean and +3 O' density values were plotted at their respective 
levels for each season of each station. A preliminary study of 
these graphs showed that the 1.86 x 10 " 2 kg m"-’ density surface 
could be expected to occur around the 30 km level. Therefore, the 
maximum and minimum height values of this density surface were 
tabulated for each season and these values were plotted according to 
the latitude of the station. Thus we get a height versus latitude 
graph showing the 99 per cent variability limits of this particular 
density surface. 

The results indicate that a small variation occurs at all 
seasons over the tropics centered at about 30 km. The T lar regions 
show a general increase in height of the surface from winter to summer 
with a very small variation during summer. The mid-latitudes show the 
same general trend as do the polar regions, but have more variation 



during summer. 

For any given season and latitude we may generalize by saying that 
the density surface will be within r 0.5 km of the mean level 99 per 
cent of the time. In fact,, considering all seasons aid all latitudes, 

—2 —Q 

we may generalize by saying that the 1.86 x 10 kg m 0 density 
surface will be within + 1 km of the 30 km surface. 

b. Height Variations of 48 km Density Surface 

In order to investigate the inde A cf ji.cal refraction 
at high altitude, a study was carried out on a strong parameter of 
this index, the density of the atmosphere. This study, in essence 
an extension of a previous study covered in a Memo dated 29 April 
I960, involves the height variation of a density surface occurring 
in the neighborhood of 150,000 feet. Due to the dearth of the data 
at this altitude, the conclusions drawn must of necessity be considered 
as preliminary and "rough". 

Various speculations have been made as to density variations at 
a given height at high altitudes but none, to the writer’s knowledge, 

g 

of the variation of the height of a given density surface. Goody 
quotes an early study by Whipple, Jacchia and Kopal to the effect 
that the upper stratosphere shows a marked seasonal variation in 
density and that the upper atmosphere does not seem to react to aay- 
tc-day changes in surface temperatures. From rocket data, Whipple^ 
concludes that there is little evidence of a diurnal '~ariation but 
that atmospheric densities are greater in summer than in winter. As 
to latitudinal variations he states that there is no indication of 
density variation between the equator and latitude 33° N but meteoric 
and acoustic data do indicate increasing densities with increasing 

h2 


latii ide. 




K.'llogg st'ves that Soviet locket observations up to 80 km show 
that the Arctic summer densities are about the same as mid-latitude 
densities; however, Artie winter densities are 10-20 per cent lower 
than mid- 1 'iude densities. Kellogg also states that, in the 
30-80 km range, tne day-to-day variations of density are greater at 
Churchill, Canada, than at White Sands, New Mexico, expecially during 
winter. 

LaGow, et al^ - studied the Arctic atmosphere using rocket data 
from Churchill (59° N)- They fina, in the 30-70 km range, that the 
summer values are 5-10 per cent higher than those derived by the Rocket 
Panel from data taken at White Sands (32° N). Also, in the 25-40 km 

range at Churchill, winter densities are lower than summer densities. 

12 

Jones, et al studied atmospheric densities using a falling 
sphere technique. Flights were made at White Sands, Wallops Island 
(38° N), Churchill and on shipboard (49-65 0 N). The White Sands and 
Wallops flights showed little r '?nsity variations. Winter densities over 
Churchill were generally lower than at White Sands. In one particular 
instance, over a 2-day period at Churchill, the 50 km density increased 
about 80 per cer.t. This was associated wita abrupt stratospheric warning. 

In summarizing the reports mentioned above, the following condi¬ 
tions seem to prevail in the upper atmosphere. The summer densities 
are greater than the winter densities. This implies the movement of an 
isopycnio xevel to higher altitudes in summer. Diurnal variations are 
slight in general, but there may be large changes in late winter or early 
spring due to abrupt stratospheric warming. Latitudinally there seems 
to be little variation in isopycnic nel^hcs from the equator up to about 
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35° N. Then the levels increase in height with increasing latitude 
up to tne polar legions where they remain constant during summer but 
show a decrease during winter. 

c. Data Calculations 

To establish more firmly the above hypothesis and 
obtain some order of magnitude of the height variations, a search 
was made for actual data from which height variations an isopvcnic 
levei could be readiiy calculated. A group of high altitude radiosonde 
flights were made over Berlin (53° N) which allowed estimates to be 
made. Values of density versus height of the 10 mb and 5 mb pressure 
surfaces were plotted on semi-log paper. For spring months there were 
42 values at the 10 mb surface and 13 values at the 5 rub surface. 

For summer months there were 40 and 19 values, and for winter there 

13 

were 27 and 9 values. Lowenthal published a group of summertime 

high-altitude radiosonde data from 7 stations. From this daca 15 

density-height values for 10 mb and higher were found for Portland, 

Oregon; 31 for Narsarssuai, Greenland, 11 for Denver, Colorado; 

30 for Belmare, New Jersey; 10 for Long Bench, California; 14 for 

Chanute, Illinois; and 11 for Goose Bay, Canada. All these values 

14 

were again plotted on semi-log paper. Finally Elterman presents 
a few values from White Sands, New Mexico (6 for summer, 12 for fall) 
obtained using a search-light technique. 

By assuming a linear relation between height and log density, 
straight line envelopes were determined for the seasons at the various 
stations. Due to the sparse data, these envelopes were assumed to con¬ 
tain plur or minus one standard deviation of density-height values. 

From these graphs the mean height and standard deviation of the 
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4.3 x 10 "° gm/enr isopycnic level was determined and plotted on a 
height versus latitude graph. This isopycnic level occurs around 
39-5 km (about 130,000 feet) and it does show an increase in height 
with increasing latitude over the mid-latitudes. Insufficient 
data for other seasons did not permit a latitudinal variation graph 
to be made. The actus? values, however, are presented below. 

Height of 4.5 x 10” gm/cnr Level 


Location 

Season 

Mean Ht (km) 

R td Bov 

White Sands 

Fall 

39-2 

0.6 

Berlin 

Spring 

38.5 

0.4 

Berlin 

Winter 

3 6.6 

0.6 


The summer standard deviations of almost all the stations 
were 0.2 km. 


The results of this study may be summarized as follows 
(keeping in mind the scarcity of data): 

1) At altitudes of the order of 130-150,000 feet, the 
isopycnic surfaces are at lower levels during the 
winter than during the summer. 

2 ) Iatitudinally the isopycnic surfaces vary little in 
the tropics, increase altitude with latitude in the 
sub-tropical and mid-latitudinal ranges and remain level 
over Polar regions during winter. .‘The increase in height 
between tropical and Polar regions appears to be the order 
of 1 - 1.5 km. 

3 ) The height variations of the isopycnic surfaces are 
smaller during summer than during the other seasons. 

During summer 99 per cent of the height values will 
fall within 2,000 feet of the mean seasonal height. 

During winter this variation increases to 6,000 feet. 

4) One reason why the spring and winter vari"vions are 
larger than the summer variation may be because of the 
abrupt stratospheric warming which occurs during late 
winter or early spring. 
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a. Short Term Fluctuations 

It is known that the equal density surfaces fluctuate 
in more or less predictable manner according to latitude and season. 

Ey correcting for the effect of these fluctuations the accuracy of 
our transit time measurement can be considerably improved. The 
question arises, however, of effect cf short-term fluctuations in 
the altitude of the isopycnic surface. 

i r- 

From meteor data it has been observed*-' that at .* mean height of 
78 km ( 256,000 feet) the total seasonal range corresponded to a height 
variation of 8.6 km ( 28,000 feet). Tf this is a more or less sinusoidally 
varying function, the rms seasonal fluctuation is about 9>000 feet. 

Ey correcting for this seasonal fluctuation, the rms error can be 
reduced to 2,000 or 3>000 feet. 

It was observed that the correlation is not improved by a com¬ 
parison with the actual ground temperature at that date rather than with 
the general average. Thus, the correlation is truly a seasonal one 
which does not measure local variations. No effects associated with 
synoptic weather fronts, deviant temperatures in the lower stratosphere, 
sunspot numbers, lunar-hour angle or solar-hour angle are conspicuous. 

From further meteor observations Jacchia^ finds evidence that the 
seasonal effect decreases with increasing height, becoming small and 
uncertain around the 100 km level. 

■The data from Jacchia supports the contention that seasonal and 
latitudinal corrections are of the greatest importance in the refine¬ 
ment of the star transit time measurements. 

Ir. comparison with these two sources of fluctuation, the random 
component or the component correlated with other physical phenomena 
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are small. In addition, our light ray traverses a considerable 
length through the atmosphere and therefore allows considerable 
averaging. For example, a ray at minimum altitude of 150,000 feet 
traverses a length of I.000 miles from the time it passes through 
the 200,Ouo foot level to zhe time it emerges from this level. 

2- Atmospheric Spectral Absorption 

In the use of spectral absorption for the definition of star 
transit times, the ozone absorption is the trosc attractive. At wave 
lengths longer than 3,000 angstroms the absorption drops off rapidly 
and the ray penetrates too low in the atmosphere before sufficient 
absorption takes place. At wave lengths shorter than 2,000 angstroms 
the 50 per cent absorption point is achieved at very high altitudes, 
which is desirable, but it involves a region of the spectrum which is 
largely unexplored and which imposes a severe restriction or. the number 
of stars usable for the measurement. 

Unfortunately, the ozone content of the atmosphere undergoes a 
wide variation from hour-to-hour, day-to-day, seascn-to-season and from 
one latitude to another. Figures (15) and (16) are illustrative of this 
point. Thus, in the case of ozone it is probably not possible to predict 
the altitude of a given isoabsorption level more accurately than 8,000 
to 10,000 feet (rms). Because the accuracy potential of determining 
altitude by refraction measurements is at least twice as good as tnis, 
and because simpler seasonal and latitudinal corrections must be made, 
ozone spectral absorption is not recommended. 

In the case of the sun, however, the use of spectral absorption is 
recommended. The sun is such an extended target that a single refrac¬ 
tion measurement cannot be made because light emanates simultaneously 
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from opposite edges of the sun's disk. Considerable accuracy is achievable 
by measuring absorption because the absorption lines of a predictable 
distributed component of the atmosphere such as COg can be used. In this 
case we will be working m the mid-infrared region of the spectrum, (in 
tht- case of tae stars, there ’"s insufficient energy in this region of Ihe 
spectrum and we are restricted to operation at wave lengths generally 
less than 1 micron.) 

A number of infrared absorption lines ior CO^ arc . iowi* m Table I. 

These lines range from 1 4 to 15 microns, and vary in usable width 

from 0.06 microns to 0.2 microns. Shown in this table is the equivalent 

sea level path for 50 per cent absorption for various wave lengths. This 

17 

data nas oeen ootainea from Altshuler. 

Because of their strength, the two most attractive absorption lines 
are those at 2-7 and 4.3 microns. The 2.7 micron line is also a point 
of high absorption of HgO - one of tne most variable components of the 
atmosphere. Because of this and because the 4.3 micron line has greater 
absorptivity, it is this latter line which is recommended. 

Figure (17) shows the magnitude of the spectral absorption at 
various altitudes for the 4.3 micron C0^ line. Fifty per cent absorption 
is obtained at a:, altitude of 160,000 feet or 49 km. At the 50 per cent 
point for radiation from a point source, an error of 1 per cent of the 
original intensity creates an equivalent error of 600 feet in altitude. 
Unfortunate 1 y, the sun is not a point source and one must resort to a 
numerical integration over an extended region of altitudes to determine 
the altitude at which its intensity is reduced to some value such as 
50 per cent. This will be discussed later. 

The data shown in Figure (17) have been derived with the aid of 
Figure (l8) assuming that the atmospheric attenuation is proportional 
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TABLE I 


INFRARED ABSORPTION LINES OF C0 2 


WAVE LENGTH 
(MICRONS) 

SEA LEVEL PATH 

FOR 50# ABS.(km) 

BANDWIDTH 

(MICRONS) 

ALTITUDE OF 
GRAZING RAY 
FOR 50$ 
ABSORPTION 1 

1.4 

3000 

+ .04 

0 

1.6 

3000 

.04 

0 

2.0 

100 

•03 

7 

2.72 

.9 

.06 

28 

4.3 

.025 

.08 

45 

4.85 

45 

.06 

10 

5-2 

200 

.06 

4 

9-32 

60 

.07 

9 

9-55 

80 

.06 

8 

10.30 

200 

.10 

4 

10.6 

250 

.05 

3 

15.0 

.0035 

.10 

54 
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SPECTRA!. ABSORPIION OF GRAZING RAY AT VAROUS ALTITIDES 
FOP WAVE LENGTH REGION, k.22 - 4-38 MICRONS 






Figure lo 



ALTITUDE - KILOMETERS 








to an exponential function of the pressure to the 3/2 power rather 
than to the air mass alone. After establishing the equivalent sea 
level path, the per cent transmission was calculated for the various 
sea level paths. 

Because the sun represents an extended source, the altitude at 
which, say, 50 per cent attenuation of its signal occurs is a function 
of the satellite altitude. This occurs because of the shape of the 
spectral absorption curve shown in Figure ( 17 ), ar.l *—a - ’ _,e of the 
effects of differential rel"action. Were it not for these effects, the 
53 per cent occultation surface would be very nearly cylindrical 
as it is for a star. The correction problem would be even more difficult 
with some ether attenuation level because the corrections would have to 
be superimposed on a conical shape. 

If atmospheric refraction is disregarded, the diameter of the sun 
at the point of tangency of the grazing ray can be written in the form 

41 . . D i, • 9.30 5 X 10-3 ^72 . 

where the apparent diameter of the sun at the earth distance is 

taken as 9 - 30 ^ x 10"3 radians 



For comparative purposes with figure (17), the values of /\ h at 
various satellite altitudes are shown in Table II. Time did not 
permit a calculation of the apparent occultation altitude as a 
function of satellite altitude; however, it is apparent that the 
correction function will be very complex for satellite altitudes 
greater than 1,000 miles. 

The sketch shown below shows the shapes of the occultation 
surfaces for various percentage transmissions f^r th s.mplified 
case where the ea^th is considered to have no atmosphere. 

50 $ cylinder 


Sun 



It was mentioned previously that the altitude at which 50 
per cent of the sclar energy (at /l = k.22 to 4.38 microns) is 
absorbed by the earth’s atmospnere varies as the altitude of the 
satellite varies. Table III shows the various correction factors 
which must be applied. 

If the sun were a point source, and if the satellite were at 
very low altitudes, the altitude at which 5° per cent attenuation is 
observed would be 160,000 feet. 


5? 





(CABLE II 

DIAMETER OF SUN AT POINT OF TANGENCY 
FOR GRAZING RAY AT VARIOUS SATELLITE ALTITUDES 


SATELLITE ALTITUDE 

S. MILES 

SOLAR DIAMETER 
RJET 

100 

44,000 

300 

77,600 

1,000 

147,300 

3,000 

257,000 

10,000 

658,000 

30,000 


100,000 

5,100,000 
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TABLE III 

ALTITUDE OF GRAZING RAY WITH 50$ SOLAR SICRGi 
TRANSMISSION FOR VARIOUS SATELLITE ALTITUDES 


SATELLITE_ CORRECTION TERMS, FEET 


ALTITUDE 

(iriles) 

Spectral 

Absorption 

Differential 

Refraction 

Refractive 

Bending 

Earth 

Curvature 

ALTITUDE FOR 
50$ INTENSITY 

100 

+700 

-rO 

-100 

-0 

160,600 

300 

+2000 

+350 

-200 

-0 

162,200 

1,000 

+4500 

+900 

-350 

>50 

165,000 

3,000 

+6800 

+3600 

-550 

-150 

169,700 

10,000 

+15800 

+6800 

-12 JO 

-850 

180,400 
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The spectral absorption correction terra can be derived from Figure 
(±7) and from Table II. The light extinction due to differential 
refraction can be derived frorr curves of the type shown in Figure (il). 
These values were obtained by taking the extinction at the center 
of the solar image and are therefore only approximate. The refractive 
bending has been derived from Figure ( 5 )• The earth curvature 
correction results from the fact that, for satellites at large 
distances from the earth, the image of the ^arin across 

the sun is curved. The last column of Table III shows the altitude 
at 50 per cent absorption including the correction terms which must 
be applied. Figure ( 19 ) shows the corrected 50 per cent absorption 
altitude as a function of satellite altitude. 

Since the uncertainty in the altitude of the isopycnic surfaces 
is about one-half mile at these altitudes, this kind of altitude 
correction need be made only for solar occultations taking place when 
the satellite altitude is greater than 500 miles. 

If possible, we would like to avoid complicating the system by 
making this altitude correction. Probably the easiest way of doing 
this is by selecting an attenuation level which is less than 50 per 
cent and whose occultation surface is (without the correction terms) 
a cone which gets smaller as the satellite altitude increases. 
Generally, this will be a linear function of the distance of the 
satellite from the point of closest approach A f the ray to the earth; 
however, it will not be a linear function of satellite a’.titude. As a 
result, the dotted curve shown in Figure (20) is obtained after assum¬ 
ing that the detector discrimination level is set for k6 per cent 
atmospheric attenuation. 

cP 

y~> 
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ELLITE ALTITUDE, MILES 





SATELLITE ALTITUDE. MILES 





Jr. deriving Figure (20) it was assumed that the shape of the cone 
mu£* be such that a convergence of 21,400 feet results when the satellite 
is at an altitude of 10,000 s. miles as shown in the sketch below- 



of the sun is 9-305 x 10 "3 radians, therefore, the discrimination level 
must be set such that 


A I = / 3-03 x IQ' 4 ■) J_ = .ohi4 
1 ( 9.305 xlO’ 3 ^ V 

Thus we will obtain a cone with bbe proper apex angle if the transmission 
level is set at about hg per cent (4 per cent less than 50 per cent). 

Under these conditions Table III m»so be modified. In tne first 
place, an attenuation level of 5^ per cent is achieved at an altitude 
cf approximately 156,000 feet. To the other correction factors (which 
will be approximately the same), we must add the conical convergence 
factor. When this is done we obtain Table IV. 

The rms deviation from a cylindrical occultation surface for a 
range of satellite altitudes out to 10,000 miles does not exceed 1,000 feet. 
By proper selection of the discriminating intensity for the solar transit 
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TABLE IV 

ALTITUDE OF GRAZING RAY WITH 
Wp SOLAR ENERGY TRANSMISSION 
FOR VARIOUS SATELLITE ALTITUDES 


SATELLITE 

CORRECTION TERMS, FEET 


ALTITUDE 

Factors from 

Conical 

Altitude for 

(miles) 

Table IV 

Convergence 

46$ Intensity 

100 

+600 

-1,400 

155,200 

300 

+2,200 

- 2,600 

155,800 

1,000 

+5,000 

-4,800 

155,800 

3,000 

+9,700 

-8,400 

157,300 

11.000 

*20,Lqo 

-21,400 

155,000 
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time measurement, ohe occulta Ij .on surface can for all practical 
purposes be made cylindrical. 
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C. Star Selection 


In the last quarterly report^, star pairs were considered which deviated 
frcm diametricity by less than 3 minutes of arc. This was a necessity for 
the rotating aberrascope since collinearizing wedges o? larger angles would 
introduce significant errors into the system. Now, for a non-rotating 
aberrascqpe, the collinearizing wedge is not required and the ascs of the 
two optical systems can deviate mechanically by the .-»»» angle as the star 
pair. The deviation from diametricity can now be as large as 2.5 degrees 
and stars of visual magnitude 3*5 or brighter can be considered. Bie 
star distance from the solar ecliptic was to be a minimum of 20 degrees. 

The use of the 3*5 visual magnitude restriction reduces the number of 
stars available frcm the general catalog to 259 from which eight pairs were 
found that satisfied the basic criteria. These star pairs are ltoted an 
pairs No. 1 through 8 in Table V. A further stipulation that the two star 
pairs be orthogonal within 20 degrees reduced the list of eight pairs to 
one pattern of four stars. 

In order to obtain other star pairs the visual magnitude restriction 
waa eased to a value of 4.7 or brighter. Using this new value a manual 
search of the Atlas Coeli 1950 . 0 , by Antonin Becvar, was made to find stars 
thi.t appeared to satisfy the prior stated conditions of diametricity and 
ecliptic distance. The total search produced a list of 35 pairs which are 
listed, to Table V. This listing includes a few pairs which deviate from 
diametricity by as large an angle as 3*3 degrees and wei-? included because 
of their desirable brightness. Of the pairs listed, those numbered down 
to pair No. 25 meet the requirement of distance from the ecliptic. 

In order to compare tho brightness of the stars under consideration, 
a correction based on temperature of the star is necessary. Ofce correction 
as a function of temperature is shown to Figure 21. A star of spectral 
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figure 21. Correction Factors for Stars 









class AO is used as the basis .or The correction; stars cf :tner 
spectral classes must have an appropriate correction added to the 
apparent magnitude in order to obtain the corrected magnitude Thus 
the magnitude corrected for type (Hp) is the sum cf the visual mag¬ 
nitude (My) and the correction factor (CT). The corrected s+ar 
magnitudes are shown in Table V along with the ether pertinent facts 
about the stars in question. 

In general, this study thoroughly covers the de and 

brighter stars, whereas those star pairs between 3-5 and 4.7 magnitude 
are studies only to a questionable extent. To determine more adequately 
star pairs that exist in magnitudes brighter than a .7 it is concluded 
that a computer study would be needed. The study would have to be 
based on stars that had the magnitude corrected for spectrum classi¬ 
fication and as such would require investigation cf stars up to a 
visual magnitude of approximately 5-0 or brighter, "itfc the magnitude 
determined and the criteria for diametricity, ecliptic distance ar.i orthc 
gonality established, an accurate list of all sta’-s meeting the desired 
conditions can be obtained. 

As a result of the present study it appears that at least six pairs 
of star" car. be found which have magnitudes equal or less than a. 3 , 
and diametricity equal or less than 3.25 degrees. These pairs are 
listed in Table VI. From this tab]t it appears that if magnitude 
is the prime consideration, the stars in group one should be used. This 
group also has pairs with deviations from diametricity of j b degrees. 
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Thus far the star se_ecticr. criteria have been such factors 
as magnituc- . spectral class, deviation from diametricity f:r the 
pair of stars, distance from the ecliptic and finally, the angle between 
two pairs of stars in the grouping of four required for the aberrascope 
system. In addition to these considerations, stars must be selected 
whose positions satisfy the geometrical requirements of the proclem. 

The basic requirement with respect to the geometry is that the stars be 
selected so as to be occulted by the earth at some ir. the 

satellite orbit. It is obvious that if no cccultaticns occur-, no 
occultaticn times can be measured. 

The necessary and sufficient condition that an cccultation of a 
star occu>- is the following: 

c 1 cos i sin £ - sin i cosiT sin - SI ) s f ( £ ) 
where c = R/a 

R = Radius of occulting sphere 
a = Radius of orbit 
i = Inclination of orbit 
S = Declination of star 
o< = Bignt ascension cf star 
SI = Angle to line cf nodes cf orbit 
The assumptions made are that the occulting body is a sphere end 
that the satellite oroit is cir alar. If we assume that the orcital 
inclination is between zero and ninety degrees and also, that the 
declination of the star in question is also between zer* and rinety degrees, 
a necessary and sufficient condition that an occuitation occur for any 
value of the angle to the line of nodes (il) is the following: 
c 2. sin (i + S ) 0 - 360° 
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STAR PAIRS . 


Mag. 

Variation 

Distance 

Correct 

I'-om 

from 

for 

Diametricity 

Ecliptic. 

Type 

( r\c 7VOOC 't 

(degrees) 

2.72 

.278 

58 

3.03 











.847 

5: 

I.O69 
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ST 1 

«i 


1- 

No. 

General 

Catalog 

No. 

R.A. 

1950 

Decl. 

1950 


12 

3300 

19777 

02 41 44.5 

14 38 45.6 

-14 04 10 

13 56 30 


13 

3484 

20226 

02 56 21 b 

15 00 3.7 

-4o 30 a 5 

40 35 12 


14 

9747 

26520 

07 16 51-6 

19 12 32.9 

-6? 51 57 

67 34 25 


15 

2289 

18845 

01 50 46.4 

13 54 oo.4 

63 25 30 

-63 26 34 


16 

2756 

19269 

02 14 43.4 

14 14 23.7 

-51 44 35 

51 35 50 


17 

3463 

20029 

02 50 46.2 

14 50 49.7 

-75 16 17 

74 21 35 

18 

4730 

21332 

03 52 59-5 

15 50 43 

62 55 4l 

-63 16 43 


19 

7287 

24221 

05 46 03.9 

17 47 52.5 

-51 05 02 

50 47 31 


20 

7587 

24509 

05 57 33-2 

17 58 08.4 

-03 04 29 

02 55 56 


21 

11959 

28780 

08 39 24.1 

20 37 18.9 

-15 45 45 

15 44 04 


22 

12923 

29848 

00 19 45.0 -62 11 28 

1 * 

21 17 23.21 62 22 24 
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PAIRS 


PABLS V 
[ cont.) 


Type 

Mag. 

Nag. 

Variation 

Distance 


Apparent 

Correct 

from 

from 



for 

Diametriclty 

r r". inti r- 



Type 

(degrees) 

(degrees) 

T>~ 

4.39 

4.29 

7 25 

80 

A2 

3.87 

3-9 1 



A2 

3-42 

3-92 

.?04 

57 

05 

3.63 

/. .23 

L 



F5 



.501 

90 

KO 





33 

3-44 

3-39 

.362 

52 

KO 

4.68 

_ 



B 8 


3.68 

•155 

65 

A5 

4.78 

4.83 



K2 

4.70 

5.42 

.912 

89 

K5 

2.24 

2.99 



R9 

4.87 

4,77 

-**35 

43 

FO 

3-04' 

3-29 



A3 

3-94 

4.04 

.400 

75 

A2 

5.09 

5.29 



KO 

4.68 

5.38 

• 152 

27 

B5p 

3-92 

3.82 



KO 

4.93 

5.68 

.506 

34 

A5 

3-36 

4.01 



KO 


5.56 

•330 

78 



2.75 
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STAR PAIRS 
TABLE V 


No. General 
Catalog 
No. 



Decl. 

1950 



Mag, ¥& g. Varr'.ati :t. Distance 

Apparent Correct from f rr,m 

for Di c 'netricity Ecliptic 
Type (degrees) (degrees) 


23 

15340 

32270 

11 06 51-6 

23 -7 32.1 

4U 46 22 

-45 31 05 

a § 

; 

3-15 

4.10 

3*85 

4.7? 

1 

.758 

24 

15547 

32450 

li 15 4f-0 

23 16 07.7 

33 22 02 

-32 48 17 

KO 

KO 

3-71 

4.51 

4.4i 

5-21 

.566 

25 

7587 

24509 

05 57 33.2 

17 58 03.4 

-03 04 29 

02 55 56 

KO 

E5 

4.68 

3.92 

5.38 

3.82 

.210 

26 

238 

16740 

00 10 39-4 

12 13 13-9 

14 54 20 

“17 15 52 

B2s 

38 

2.87 

2.78 

2.83 

2.68 

2 .434 

27 

2538 

19029 

02 04 20.9 

14 03 31-1 

23 13 37 

-26 26 33 

gK 2 

gK3 

2.23 

3-48 

r .98 

1.26 

3.220 

28 

6029 

22640 

04 53 44.0 

16 46 55-2 

33 05 20 

*■34 12 16 

gK3 

gG9 

2.90 

2.38 

3.65 

3 06 

1 . 8 -J- 

29 

8208 

25180 

06 19 56.1 

18 24 53-0 

22 32 23 

-25 2? 04 

gM3 

gKl 

3-19 

2.94 

^ <01 

3 

3.120 

30 

3391 

19954 

02 47 02.1 

14 47 20.7 

27 03 20 

-27 45 12 

B8 

K2 

3.68 

4.63 

3.58 

5-38 

.700 

31 

5599 

22332 

04 32 54.3 

16 34 24.1 

10 03 35 

-10 28 03 

A3 

BO 

4 *33 

2,70 

4.48 

2.65 , 

.549 

32 

6306 

23158 

05 06 50.1 

17 07 30.5 

15 32 06 

-15 39 53 

FO 

A2 

4.86 

2.63 

5-10 

2.68 

.210 
















































































































STAR PAIRS 


TABLE V 
(cont.) 


No. 

General 

Catalog 

No. 

R.A. 

1950 

Decl. 

1950 

Type 

M3g- 

Apparent 

Mag. 

Corireo; 

for 

Type 

Variation 

from 

Diametricity 

(degrees) 

Distance 

from 

Eccliptic 
(degrees) 


155.11 

11 14 07.1 

-03 22 4l 

A5 


4.73 



33 







.386 

8 


32415 

23 14 34-3 

03 00 32 

mm 


4.55 




16425 

11 58 18.6 

06 53 35 

A3 

4.57 




34 







.658 

7 


33330 

23 59 23.7 

06 17 31 

M6 

4.66 





163A9 

11 46 30.6 

14 51 06 

A2 

2.23 

2.27 



35 







1.540 

12 


32931 

23 40 07.8 

-14 49 18 

AO 

4.62 

4.62 
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TABLE VI 


AEERRASCOPE STAR GROUPS 


Group 

General 

Constellation 

Type 

Magnitude 

Bwi&bicr. 

Angie 

Distance 

No. 

Catalog 



(corrected) 

from 

between 

from 


No. 




Diametricity 

Paira 

Ecliptic 






(degrees) 

(degrees) 

(degrees) 

■ 

2572 

j3 Tri 

A5 

3.09 

1.5 


23 


19033 

9 Cen 

KO 

2.66 




H 

53.64 

ot Ret 

G5 

4.06 

1.5 

102 

84 


22101 

Dra 

G5 

3-59 





645 

Caa 

cBOe 

4.20 

3-2 


59 


17374 

fi Cru 

B1 

1.46 




2 

3300 

nf Cet 

B5 

4.29 

0.7 

99 

30 


19777 

J Boo 

A2 

3-91 





l4oo 

fit And 

MO 

3-22 




3 

18039 

i Cen 

A2 

2.97 

2.4 


28 

3584 

9 Eri 

A2 

3-92 

0.7 

99 

57 



20226 

B Boo 

G5 

4.23 
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For values of the angle to the line of nodes between the right ascen¬ 
sion of the star and the right ascension plus 45°. or in the range 
right ascension plus 135 ° to right ascension plus 360 °, the following 
is a necessary and sufficient condition that occultation occur: 

c ^ cos i sin £ + ( 1 / \[2 ) sin i cos S 

+ 45 ° 

od -r 135°-'.°, < '~>i * 360° 

Thus for 75$ of all values of .0. , an occultation will occur. 

For values of the angle to the line of nodes between the right 
ascension pins l 80 c and right nsr^nsion plus 360 °, the following is 
a necessary and sufficient condition that occultation occur: 

c rh cos i sin £ 

In this case, for 50$ of all values of Si , an occultation will occur. 

These three conditions are illustrated graphically in Figures 22, 

23j and 2k As an example, (for a 3 00 mile orbit where c = .93) 
assume that an orbital inclination of 45° is desired. Then, if 
occultation is necessary for all values of Si. , it can be seen that 
stars must be selected which have a declination between 0 and 23 " 
or between 67 " and 90 ° where .93 > f (<£). 

The reason for considering varying values of the angle to the 
line of nodes (-ft) is the regression of the line of nodes with time. 

The >*ate of regression (Ait) is given by the following:^ 

AfL= —£££ A radians/rev. 

2 d 

a(l-e ) 

where J is a constant (1.633 x 10*^) and "e" is the eccentricity of the 
orbit. This equation can be used to determine the time in years required 
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Declination of Star 
Figure 22 

STAR DECUNAT1CW SELECTION FOR OCCUITATICN FOB ANY ANGLE 
TO THE LINE OF NODES OF THE CEBIT (A) 
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Declination of Star 


Figure 23 

STAR DECLINATION SELECTION FOR OCCULTATION FOR ANGLE 
TO TEE LINE OF NODES (/2.) AS FOLLOWS: 

e*+ <joc +360’ 
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for the line of nodes to regress thru an angle of 360*. For a circular 
orbit where "e" equals zero, this tine in years (T) is expressed as 
follows: 

T = a 5/ 2 secant i x 10”9 

This equation is plotted in Figure 25* Thus for a 300 mile orbit at 
zero inclination, 1.3 years will be required for thr line of nodes to 
regress 360*. 




Radius o£ Circular Orbit (a) miles 
Figure 25 

TIME FOR REGRESSION VS. RADIUS OF ORBIT 










IV SYSTEM ERROR .ANALYSIS 


In this section two related and complimentary topics are discussed. 
Hie first concerns those erors introduced into the final calculation by 
the effects of meteorological irregularities and instrument errors; and 
the second concerns computational errors. Errors of the first topic 
result in errors in the measured time of occultatiOn. .‘Ithough these 
errors can be known and compensated, errors will still remain in the 
final calculation because of the facts that the system equation does 
not perfectly fit the problem, and the computer imposes limitation on 
the accuracy with which the problem can be solved. 

A. Physical Errors 

The sources of errors can be enumerated as follows: 

1) Meteorological Errors 

a) Systematic errors in mean latitudinal and seasonal 
corrections to nominal isorefraction altitude. 

b) Statistical fluctuations in atmosphere causing 
deviations from assumed altitude. 

c) Scintillation due to small scale inhomogeneities 
in the atmosphere resulting from wind shear. 

d) Background noise due to scattered radiation from 
the sun. 

2) Instrument Errors 

a) Systematic errors in measurement of the time at which 
refraction builds up to a preassigned value. (This 
error is equivalent to an angle error in the instru¬ 
ment refraction measurement). 

b) Random instrument errors in transit time measurement. 

c) Error due to changing aberration during refraction 
build-up. 
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d) Error due to altitude drift of platform because 
stars are not perfectly diametrically opposite. 

e) Error due to reference clock. 

3) Computational Errors 

Given a sufficiently large computer, the computational 
errors should be much smaller than those listed under (l) and (2); 
however, since our goal is a self-contained positional accuracy of 
about 1 mile achievable with a computer of minimum we'--'.-. certain 
sources of computational errors will require investigation. We 
may not be able to ignore the ellipsoidal shape of the earth, the 
nodal regression, the precession of the perigee, and the problem of the 
rate of convergence of the successive iterations. 

To indicate the influence of some of these errors, it ■>& instruc¬ 
tive to consider the sensitivity of the satellite altitude claculation 
to the errors in the altitude of the occultation surface. To do this 


we will examine a simplified case in which the plane of the orbit is 
known and in which the orbit is circular. If we select two diametrically 
opposed stars which lie in the orbital plane, a geometry of the type 
shown in Figure 26 results. (While this is valid here, see conclusions 
of paragraph IV B concerning in-plane stars.) The parameters shown in 
Figure 26 are listed below: 


j*" = Orbital radius 

J> 0 = Radius of occulting surface 

J* e = Radius of earth 


= Earth's central angle which subtends one half 
of the arc between ingress and egress 
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I 


Fig-ore 26 

Xc-Plare Star Fairs 
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If t = time of egress of star 2, and 
e 

t. = time of ingress of star 1, then 


Q 

ir ’ 


therefore. 


.■co s tLlS : !r > = £° 


T is a function of P but P is an implicit function of P , t , and 
J s -• s J o e 

t . Thus, we cannot obtain a direct solution for P . However, we can 
i 's’ 

obtain the error sensitivity coefficients in explicit form. 

If we assume that the errors in P , t , and t. are small and un- 

J o’ e i 

correlated, then we need retain only the first order term in tne Taylor's 
expansion and may write 

<r iff - iff * tpfr + • 

This amounts to separating the meteorological and instrument errors, i.e., 
the (t) are rms errors due only to the instrument time measurement and 
are not functions of the isorefraction surface errors, O'{f Q ) • To 
determine the value of we take the partial derivative of the cos )( 

3/\> 

equation with respect toy^. Thus 


'Tf T . y 1 — s ) 

- 2 — sin # = t ' 

T ofo ys J s dfo 


From the expression for the orbital period we can write 


d_i_ = jl 

'dfo d/’s ofo 


8i 







mwau 


Therefore 


A 


J /?( 


4 s +1 * sin % 

r s 



In the limiting cases we have lim 

Vs 

< A 

1 

f 3 b"Vo 

1 

■c 

lim 

(A 

or/^O 

] 

\ 

/A"* ^ 


If 

3tt 


.h 25 


3 /? 

Because of the maxima value = 1, we can conclude that the effect 

a /0 o 

of errors in the assumed altitude of the occultation surface is to 
create errors in the calculated satellite altitude which are no greater 
than the f Q errors. Actually, If we had completed one orbital rotation 
and measured T directly would not be a function of /° and therefore 

<r J ^ fi 

- r - ■~ a * 0* Figure 27 shows the manner in which —_ZJL varies with /?. 
Vo 

If venow take the partial derivative with respect to t e (or tj_) in 
the same manner as before we obtain the expressions 

[t - (t- t*) (-y~)7 ®in y ■ -Vfg , ( j A -) ; 

L * Tte-1 7T/>s ^ L t e 


^ T 


3 T / f b \ 
Vo ( ^e 


Therefore 


Vs 




[ 


sin 


■4 s + | ^sin ^ 
/ s 2 


] 
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where 



The value of this error sensitivity coefficient under the limiting 
cases hecomes 



11m 

/s" 6 * 



0 


An examination of the original equation will show that 



Figure 28 shows the values of 



for various satellite altitudes. 


Note that this sensitivity coefficient has a maximum. Bie location of 


this maximum has not "been evaluated analytically, 
of Figure 28 will reveal that the maximum value of 


however, an inspection 


Jjfe 

<)te 


.96, and that 


this maximum lies at a satellite altitude of about 1200 miles. 


One or two comments are in order regarding the implications of 
Figures 27 and 28 on the ultimate accuracy achieved with a system of 
this type. Since the rms errors in /° from meteorological uncertainties 
are about 0.5 mile, Figure 27 indicates that the errors in the satellite 
altitude due to this effect alone will be equal to or less than 0.5 mile. 
If we allow the refraction to build up to 10 seconds of arc, then the 
refraction rate will range from 2 to 6 seconds of arc per second. It 
is not unreasonable to expect that the refraction can be read to an 
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ERROR SENSITIVITY COEFFICIENT, . .. ■ 8 AT VARIOUS SATELLITE ALT3 FUDES 
















accuracy of at least 0.5 seconds of arc and therefore the instrument time 
error will not exceed .08 - .25 seconds. Under these conditions the 
maximum error An due to instrument errors alone will be about .2 
miles. 

It therefore appears that nhe errors due to the meteorological 
fluctuations will be larger than the errors due the instrument re¬ 
fraction errors. By determining the semi-major axis from one complete 
orbital rotation the effect of some of these meteorological uncertainties 
can be reduced. If the isorefraction surface follows some kind of reason¬ 
able function (of perhaps 2 or 3 parameters) and if the statistical 
fluctuations from this mathematical mode is small, then additional star 
transits can be used to effect a large improvement in the accuracy of the 
system compared with the case where one must estimate the shape of the 
iaoraf.ractiou surface from prior data. 

In some applications it is necessary to be able to calculate the 
satellite velocity with great accuracy. Since 



If/^ k 396 O miles, b a 3960 + 300 miles, g * 32 ft/sec, then 
the error sensitivity coefficient is 

■ = 5 - 5 ^ x 10 "^ feet/sec/foot 

As was shown above, the errors in g will be of the order of one mile; 
therefore, the error in the calculated satellite velocity will be of the 
order of 3 feet/sec. 
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B. SYSTEM EQUATIONS FOR THE OCCUHTATION TECHNIQUE 


The derivation of the equations upon vhich the occulta+don 
technique is based will be sketched in this paragraph. Two 
assumptions are first made which, while considerably simplifying the 
analysis, do not seriously limit the applicability of the technique. 

These are: 

1) The occulting body or planet is considered a sphere. 

Thus, the shadow boundary formed by the stars and the 
planet will be a right circular cylinder extending 
from the edges of the planet to infinity in a direction 
away from the star. 

2) The second assumption is dynamic and states that the 
orbit shall be a Kepler eclipse in inertial space, i.e., 
the satellite moves in a central force field. 

The assumption (2) can be relaxed somewhat to allowing effects 
such as nodal regression, apsidal precession, etc., which are sufficiently 
small. This is valid since the proposed technique gathers sufficient 
data for a complete orbital determination in one satellite revolution. 
Although the mathematics become somewhat more complicated, assumption 
(l) can also be amended to include planets in the form of oblate spheroids; 
however, this refinement is not discussed here. 

We now proceed with the analysis. Consider a cartesian coordinate 
(X, Y, Z) system with origin at the center of the earth, Z extending along 
the north pole, and X pointing toward the first point of .'Vries. This 
is called the "astronomical frame," or the inertial frame, and it is the 
one with respect to which the parameters of satellites are usually given. 
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If we employ another reference (x, y, z) frame having the same origin, 
hut with the x-y plane containing the plane of the satellite orbit, we 
can write the position of the satellite with respect to the astronomical 
frame as: 


where 


and 


with 



/ B 



\ 

r 

* 

/ 


_ a (1 - e 2 ) 

/ - 1 + e cos 9 

= instantaneous radiuB of the satellite 
a c semi-major axis of the orbit eclipse 
e a eccentricity 
CJ - argument of perigee 
0 •= real anomoly (function of time) 


( cos XL. - cob i sinsin I sin -fl \ 

siu/L cos i cos jx ain i cos A- ) 

0 sin i cos i J 

-fL = argument of the line cf nodes, and 
i = inclination of the orbit plane 


Furthermore, under assumption (1), we can construct a third 
cartesian reference system (x£, y£, z£) again having the seme origin 
but with extending toward the star under discussion (the k^ star), 
and y£ lying in the X-Y plane. In this frame the equation of the 
occulting right circular cylinder becomes 
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where 



/Vo C0 Vk \ 

! /o 3i Vk 


\ 


/ 


/0 Q si radius of the occulting body, and 

a angle in the x£ - y£ plane measured from the x£ axis 


Not at the moment of occultation, some point on the orbit, eclipse is 
coincident with some point on the occulting cylinder. To mathematically 
obtain the implications of this physical fact, we must transform (l) 
into the k th stellar frame (x£, y£, z£) and set it equal to (2). Hie 
transformation is accomplished by means of the orthogonal matrix: 


/cosa k 

sino6 k 

cos & k sin d k 

- cos<J k 

-sinoi k 


cosol k 

0 

cosol k 

cos^ k 

8inc< k cos £" k 

sin £ k 

right ascension of the k"^ star, and 



where 

*k 

£ k ■ declination of same 

Thus, the final geometrical relation is embodied inthe matrix equation: 


Vo *»M i\ 


A> 8i Vkj 


/ 



(3) 
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Here the subscript j can take on values 1 and 2 corresponding to ingress 
and egress respectively. 

Before obtaining the final explicit fora of the equations for the 
occultation technique, let us examine equation (3) in the ligit of what 
can be experimentally measured. It is currently proposed to measure two 
quantities during an occultation. One of these is the time at which the 
magnitude of the refraction becomes 10 sec while the other in the direction 
of travel of the ray in inertial space as determined >>j the ratio of the 
two components of refraction. The latter is measured in a plane perpendi¬ 
cular to the sight line to the k"** star and therefore constitutes a 
determination The former can be ’iaed as an input to €».. Thus, 

we can only vuso two of than and we shall now write these in a form more 
amenable to the inputs. Upon squaring and adding all three equations of 
(3) we find that 

A 2 - A 2 ♦ <“) 

! This yields a relation independent of k j : 

\ A 2 r f 

1 -"7^2 r [ COB h 

| + COB 6 k 

r +sin£ k 

i 

An equation independent of both f q and ^ can be obtained by dividing 
the second equation of (3) by the first: 


cos (e. +y) cos (a, - xi) 

J £ 

sin (8j +cj) cos i sin (<* k - - 4 l) 

sin i sin (Oj +W )J 2 (5) 
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= £cos i sin (@j + 10 ) cos (CX^ - A.) 

-cos (Gj+CJ) sin (« k 
X jjsin cos (©j + u>) cos (oj^ -Jh) 

+ sin <£p_ sin (©j + oJ ) cos i sin (<X k - JX) 

-cos £ k sin i sin (©^ + oJ )~j (6) 


These latter two equations form the basis for the occultation technique. 

The known quantities are usually <x k> £sod tj^: the time 

of occultation. The unknowns are a, e, i, w, , and T 0 : the time at 
perigee (contained on ©j). Since we obtain one pair of equations (5) and 
(6) at each occultation, it is easily seen that 3 occultations are needed 
to completely determine the 6 orbital parameters. 

In the case of solar occultations, it does not seem feasible to 
determine M kj and therefore in this case, equation (5) would be the only 
one usable, in this instance, one then only needs 2 stellar occultations 
and 2 solar occultations to determine the parameters of the satellite orbit. 
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C. Results of Error Analysis 

An error analysis was performed during the last few months 
on the equations basic to the occultation technique. In order to do this, 
specific orbits and stars had to be chosen, but the results seem to lend 
themselves to general conclusions regarding the applicability of the 
technique. 

Hie stars chosen, together with some associated characteristics 
are listed in paragraph IIC of the first quarterly report. These form 
a tetrad consisting of 2 pairs of diametrically opposed stars approximately 
orthogonal to each other. Hie inclination which their plane makes with the 
equator is about 60 °. 

There were several sets of orbits chosen, but those which were 
examined closest are discussed here and exemplify the general results. 

Hie class of orbits examined herein all have the following parameter in 
common: 

a = 4260 miles 
e = .01 
i = /2 

OJ =0 

T = 0 
o 

The omitted parameter fi. was allowed to take on 6 values distributed at 
60° intervals about a circle. 

While the full technique proposes to measure the time as well as 
the direction of the stellar ray during occultation, this analysis assumes 
that only time is measured. Furthermore, with the class of orbits considered, 
there was no jTl for which one of the stars did cot occult. Thus, during 
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each rotation, 8 times -were measured. This represents an overdetermination 
of the results. Since only 6 inputs are needed to calculate all of the 
orbit parameters there are 8.'/ (8 - 6) .' 6 ! =28 possible combinations 

from which to do this, lb further complicate matters, it was decided to 
calculate the semi-ma/jor axis a from the period measurement: 

£ - ^ a 3 
k 

) 

Ihen the error in a amounts only to 0.47 miles for an error in the period 
of one second of time. Since this is considerably better than is expected 
from the full occultation technique, we will assume that a is given in 
this manner. Then only 5 parameters remain giving 8! / (8 - 5) •' 5 1=56 
possible combinations to use. Therefore a rather arbitrary selection of 
inputs was made according to Table VII. It might be noted here that it . 
will be essential to select the best numerical data in the presence of 
redundant but independent information (if the data is dependent, the selection 
is not so difficult as discussed later). While a weighted average is often 
helpful, this will probably not be as useful in the present problem as 
the development of criteria for the selection of those measurements which 
can yield the least rms error in the final calculation of the satellite 
position. As a result of our instrument error analysis the most accurate 
measurements will be as follows in order of accuracy. 

a) star transit times, 1-2 miles 

b) solar transit times, 2-4 miles, 

c) star refraction rate, 10-20 miles 

d) star refraction direction, 20-40 miles, 

e) star intensity 30-60 miles 
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TAmLE VII 
INPUT SELECTION 


7- 

3T/iK 

INGRESS 

EGRESS 

1 

X 

X 

2 

x 

X 

-a 

X 


k 




The calculation proceeded as fellows: the proc^r *n the AD/E-CS* 
\ F 

computed all the —2^ - 


where x. are the orbit parameters. Then the inverse of the resulting 

* F k 


5x5 matrix multiplied by a diagonal matrix consisting of 


produced all 36 of the sensitivity coefficients 




t. 


A 


6 ^ 


where are the occultation time measurements. What is wanted in the 


final analysis is the sigma error; 

2 \ 2 
Xl 

i» l 




(x^) are plotted ir. Figures 25 through 33 against 12. assuming that 
o'MJ =cr(t 2 ) = . . . = cr(t 5 ) ^ .1 sec. 


Of course, the plot of o'M vs -SI does not appear since this is a 
•^nstart due to its method of derivation from the orbital period. 

These graphs have at least one thing in common; two peaks situated at 
50° and 240°. Upon examining the geometry it is found that in both 
of these instances, the line between the star pair 1 and 2 lies nearly 
in the plane of the satellite orbit. Since, as seen from Table VII nearly 


^General Mills Computer, Adaptable Digital Electronic Computer System 





cr(s) vs. J1 
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all of the inputs come from this pair, the general conclusion can be drawn 
that stars which lie in the plane of the orbit lead to large errors in 
the determination of that plane. Biis conclusion is supported by the 
other classes of orbits considered, though not discussed here. 

Another conclusion can be obtained by considering the relationship 
between the eccentricity of the orbit and the diametricity of the stars. 

If ve select a circular orbit, i.e. one with e = 0, the-i., since the 
occulting surface cuts the orbit plane in an elipse, ve have the relation 
shown in the figure below: 



It can easily be seen that a measurement of the times at 1 and 2 
leads uniquely to a determination of the times at position 3 and 4. 

Bius the equation derived for these are dependent upon those for position 
1 and 2 and hence do not yield additional data, but merely repetitive data. 
Ihis is a consequence of the diametricity of the star pair coupled with 
the circular shape of the orbit. To apply this to our case, we see 
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that the orbital eccentricity e = .01 is quite small yielding a nearly 
circular orbit -while the stars are within a few minutes of arc of being 
diametrically opposite, Hius from the above discussion, an explanation 
of the large errors encountered in the graphs is attained. Note that 
this does not only apply for the two peaks ati)_= 60 ° and 2hO°. Hiis 
is a general effect dependent only on the circularity of the orbit and 
diametricity of the star pairs. 

In the light of these observations, three other points were calculated 
and plotted on the graphs. In all of these it was assumed that a, e, 
and ui are given from other sources in order to do the calculations by 
hand. In other words, only cr (i), CT (/2), and <r (t ) were found, 
assuming the remaining to be zero. Also the interaction between the 
low eccentricity of the orbit and the near diametricity of the star 
pairs was taken into account in the selection of input data. 

Hie point at-A - 2kO° denoted by [[} used inputs: 

Ingress Egress 

X X 

X 

It is seen that the errors are indeed smaller, but still do not become 
as low as desired. However, this is the very value offlfor which star 
pair 1 and 2 lie in the orbit plane. Ihus lhat should be investigated 
is an/lvalue not 60° or 240°. 
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The point at /L= 300° denoted by £TJ satisfies this somewhat. The 
same inputs were used and the errors also are reduced to some extent; 
however, while pair 1 and 2 are not in the orbit plane, pair 3 and 4 are 
fairly close and hence may contribute to the error. 

The point at Y). = 240° labeled ^ used inputs: 

Star IngresB Egress 

1 
2 

3 X X 

4 

This star does not lie in this plane of the orbit, nor do we use 
dependent information from diametrically opposite stars. Only two 
pieces of data are therefore available at this value of Y2. . Therefore, 
it was decided to calculate only Yl. and i using T Q from some other 
source. It is seen that in this instance, (Y2.) and < 7 - (i) are 
considerably reduced though still not as far a.' we would like in view 
of the possibility of <r~ (t) being greater than .1 second. However, 
it must be noted that star # 3 > while not in the orbit plane, is not dis¬ 
placed considerably from it. The angular distance is about 33 0 . 

In summary, the conclusions which can be drawn from the analysis are 
as follows: 

1) It seems most important that information derived from 
occultations of stars lying in the orbit plane be not 
used to calculate this plane. 

2) If the satellite orbit is nearly circular, information 
from diametrically opposed star pairs, even though they 
be out of the orbital plane is redundant and cannot 
therefore be used to generate a set of simultaneous 
equations which are uniquely solvable. 
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In spite of these restrictions, if we have at least one 
star pair out of the plane of the orbit, and we make use of solar 
occultations in the necessary eases, we can determine with the 
required accuracy all of the orbital parameters. For instance, 
if one star pair line lies in the orbit plane and one stellar egress 
occurs in a sun-lit atmosphere (thereby introducing sufficient 
noise to render this time measurement unusable), the semi-major 
axis, a, can be determined by a period measurement using tne in¬ 
plane stars while the remaining 5 parameters can be calculated from 
the 3 remaining out-of-plane stellar occultations plus two solar 
occultations. If the orbit is circular, we have more occultations 
than unknown parameters and can therefore select the combination 
that produces the most accurate results. 
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V. INSTRUMENT DESIGN 


A. General Description 

In the previous quarterly report, emphasis was placed on the 
analysis of an instrument which was capable of measuring aberration 
angles, and from these measurements, the position of the satellite 
over the earth was computed At the same time, consideration was 
given to the phenomena of refraction and stellar occultaticn. It 
was noted that by measuring the time at which refract-on built up to 
a prescribed level and by ascribing this time as the time of occultation, 
a very accurate determination of the orbital period can be made. Ey 
noting occultations of several stars in sequence, the orbital elements 
can be determined by the solution of five simultaneous equations 
with five unknowns. 

Since the measurement of the refraction phenomena and the subsequent 
occultation time holds promise of more accurate position determination 
than the aberration measurement, GMI proposes to build a satellite 
guidance instrument based on the refraction measuring capability. The 
proposed instrument will be capable of measuring aserration as well 
as refraction; however, since the refraction rate is much greater 
than the aberration rate, the instrument will have a higher response 
rate but will require ^ess accuracy. An occultation time measuring 
instrument will require to measure the time at which the angle 
between two stars changes to some arbitrary value. As the line of 
sight to a particular star dips down into the atmosphere of the earth, 
this sight line will be refracted to a magnitude of 10 seconds of arc 
in a matter of only 10 seconds of time. The arbitrarily selected value 
will fail within this range of refraction. , 
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The star tracker arrangement which is proposed for the occulta- 
tion technique is similar to that previously proposed for the aberra¬ 
tion technique in the following respects: 

1) The complete system is composed of 4 star trackers. 

2) Two pairs of diametrically opposed stars are used. 

3) The star pairs are approximately orthogonal to one 
another. 

h) The measurement of the relative positio--i of the 
stars is independent of the direction ci pointing 
of the telescopes. 

5 ) Photomultipliers are used for radiation detection. 

6) Two Cassigrainian telescopes are mounted back-to-back. 

7 ) A knife edge prism to be used for image splitting. 

The star tracker arrangement which is proposed for the occultation 
technique differs from that previously proposed for the aberration 
technique in the following respects: 

1) No rotation of the aberrascope is used. 

2) No optical elements external to the telescopes are 

used (such as collinearizing wedges or Herschel wedges). 

3 ) The diametricity requirement is ~elaxcd. 

» 

4) Brighter stars can be found because of the relaxation 
of the diametricity requirement. 

•» 

5) The ultimate accuracy requirement in the measurement 
of relative star positions is relaxed by one order of 
magnitude. 

6) The response rate of the error detector is increased 
by approximately one order of magnitude. 

7 ) A smaller, lighter weight optical system can be used. 

8) The electronics can be simplified. 

It is recognized that under certain restricted geometries, the 
occultation technique does not supply adequate information. An 
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investigation of other sources of information has indicated that these 
problems can be avoided. Thus the occultation technique provides almost 
as general a solution for the determination of the orbital elements as 
the aberration technique> 

The important concept concerning the occultation measuring instru¬ 
ment is the fact that instrument stability is required over a very short 
period of time, a matter of 15 seconds at best-. This means that the 
absolute alignment of the star trackers relative to one another is less 
important and that the effects of thermal transients are small. When 
refraction measurements are not being made the instrument is used for 
attitude control which can be maintained to an accuracy of at least 
5 seconds of arc The absolute alignment errors of the optical system 
(between the opposed telescopes) are not significant and rotation of the 
system to cancel out these errors will not be necessary. This change 
from the aberration measuring instrument is very advantageous in that 
the problem of designing a large bearing of high, accuracy to operate 
in the hard vacuum environmen' is eliminated. Also, the requirement 
for slip rings to provide power and taKe out electrical signals is 
removed With these changes, it is easy to see that a substantial 
simplification has been effected in the instrument design. 

Now since aberration is not being measured, the requirement for 
the Herschel Wedges is no longer present. An angle measurement is 
still necessary, but a higher response rate will be necessary. The 
Herschel Wedges, b^ing massive devices, could not conceivably be 
considered for this use. Extremely high drive forces would be neces¬ 
sary to accelerate the wedges and the resultant dynamic error would 
be large Without the Herschel Weoges, the burden of making angle 
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measurements falls upon measuring the displacement of the center 
of the knife edge detector. Tftis alternate van mentioned in our 
Quarterly Report^ ana represents a clean design, provided that 
adequate accuracy can he preserved. The use of an interferometer 
is one possible method by which this measurement can be made. 

Since the entire instrument will not be rotated, the require¬ 
ment for diametricity can be relaxed. Without the rotational feature, 
the knife edge detector is not required to be driven ;hru an angle 
proportional to the deviation from diametricity with each revolution. 
The collinearizing wedge can also be disposed and stars can be 
selected which deviate from diametricity by larger angles, for 
example, 2 degrees. This in turn will permit the selection of 
brighter stars and the aperture of the optical system will be 
reduced to 4 inches. It will also be shortened from about 32 
Inches to 24 inches in over-all length. 

A single axis detector (a knife edge prism) will be used for 
the error detector of the instrument. Since the entire system is 
non-rotating, it will now be necessary to rotate the knife edge 
detector. This can be done quite readily by building the knife 
edges on the ends of the shaft of a small motor. This motor will 
be hermetically sealed with its own atmosphere and windows will 
be provided on each end of the package to admit the incoming stellar 
radiation. The rotor of the motor vil.1 be mounted on air bearings so 
as to remove the possibility of oil settling on the optical surfaces. 
This rotating error detector will be the heart of the refraction 
measuring instrument. 
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The stellar radiation which impinges on the rotating error 
detector will he reflected outward in a conical pattern. This energy 
must he collected and refocused on a photomultiplier tube. The use 
of a section of an ellipsoid of revolution will accomplish the 
necessary collection of energy and chopping can he performed by 
making alternate portions of the ellipsoid black or mirrored surfaces. 

Thus the image is rotating and the chopper is stationary; a reversal 
of the usual procedure but nonetheless an adequate ^iutlon. 

The signs 13 from the photonr '.tiplj. ?rs will be combined in 
suitable logic circuitry and appropriate error signals will be 
produced. These error signals will be used to drive four servo motors 
which will articulate the ends of th? rotating error detector package in 
directions normal to the optical axis. The servo motors (two on each end, 
one for up-down motion and one for sidewise motion) will position the 
rotating error detector so as reach a null position where the in¬ 
coming image is centered on the knife edge. The magnitude of the 
linear motions required at the ends of the rotating knife edge package 
is plus or minus 0.050 Inches. This figure is obtained by assuming 
an effective focal length (EFL) for the optical system of 100 inches; 
then allowances for the following motions are made: 

Annual aberration 20 seconds 

Satellite aberration 10 seconds 

Attitude errors 10 seconds 

Refraction 10 seconde 

Total 50 seconds 

To be conservative and also to allow for acquisition, the maximum 
travel will be assumed to be more than twice the value totaled above 
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Allowing for both positive and negative values, a total range of 400 
seconds of arc can be achieved which in turn is '•quivalent +o a 
linear motion of 0.20 inches or plus or minus 0.10 inches. 

The servo drive mechanisms which will articulate the rotating 
error detector package will be hermetically sealed units. Each unit 
will contain a servo motor, tachometer, gearhead, and mechanism for 
converting the rotary output into a linear output. With the linear 
motion requirement o:uy 0.20 inches, the output motion can be tra- 
mitted thru a bellow or a diaphragm thus permitting the rotary com¬ 
ponents to be sealed from the outer environment. Using th; size 8 
servo components, the entire sealed unit need not be la-,.er than 
one inch in diameter and four inches in length. 

The output ij-om the non-rotating aberrascope will come from 
two en'.onew.s which will measure the translation of the center of the 
rotating error detector package. These encoders could conceivably 
oe interferometers as already mentioned where displacements of 
between 5 and 10 millionths of an inch can Le sensed. This can be 
done by using the blue or violet lines of helium light, noting of 
course, that interference produces dark bands wnerever the distance 
between the reflecting surfaces is one-half wave length or a multiple 
thereof. The techniques of interferometry are known arts and we feel 
that this approach can be implemented in the non-rotating aberrascope 
i" a straight forward manner to produce the desired outputs. 

From the interferometer outputs, some digital computing type 
logic networks will be necessary tc finally determine the required 
parameters which win produce the occultation time as required. 
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B< Optical System 

1. Structural Analysis 

In previous reports the mechanical structural requirements of 

the aberrascope and possible sources of aberrascope structural distur- 

2k 

bance were considered. The present discussion is an analysis of an 
assumed structure for the rotating portion of the aberrascope. The 
structural requirements for a rotating aberrascope are much more 
stringent than the requirements for a non-rotating instrument. Thus, 
if the requirements of a rotating instrument can be met, no proble n 
nil exist with the design of a new non-rotating instrument. 

The mechanical structural analysis has two objectives: First, to 
determine feasibility of designing a structure of practical weight which 
has sufficient rigidity to meet the accuracy requirement; and second, to 
establish reasonable estimates of structural design parameters to serve as 
starting points for more detailed analysis during actual instrument design. 

In previous discussions the sources of aberrascope structural dis¬ 
turbance were classified according to their origin in the aberrascope, the 
satellite, or the orbital environment and according to the manner in which 
the disturbing forces would be generated. In connection with the structural 
analysis these disturbing forces may be grouped in a more general way into 
the following classification: 

1) Static forces 

2 ) Vibration 

3) Thermal effects 
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Static forces may be either actual static loads or forces which vary 
slowly relative tc "uiiG lowest resonant frequency of the structure. Examples 
of such slowly varying forces are acceleration forces due to unbalance 
of the rotating aberrascope about its axis of rotation, bearing eccentricity 
and wobble, and nonorthogonality of bearing surface and aberrascope 
geometric axis. Acceleration forces from these sources will be periodic 
with the rotation of the aberrascope structure. Since the speed of ro¬ 
tation will be very low compared to the lowest resonant frequency of the 
structure, the structural deflection under these slowly varying loads will 
be essentially the same as under static loads of an equal magnitude^ 

Dynamic loading of the structure may result from vibrations. Examples 
of possible vibration sources are bearing and gear mechanical noisa, drive 
motor vibrations, and vibrations due to moving parts or rotating elements 
in the satellite. If the frequency of the disturbing vibration is near a 
resonant frequency of the structure, a large dynamic force amplification 
will occur. In a lightly damped structure, such as the aberrascope, 
amplification factors in the range of 25 to 50 may be obtained at resonance. 

Deflections can occur because of thermal effects on the aberrascope 
structure. The heat sources may be in the aberrascope, in the satellite 
or in the orbit environment. 

The aberrascope structure may be considered to be divided into an 
inner structure and an outer structure. The inner structure is the ro¬ 
tating part of -file aberrascope which supports the primary mirror, secondary 
mirror, and the knife edge bar. The outer aberrascope structure provides 
attachment to the satellite and interconnection between crossed aber.ra- 
scopes. 
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The structural analysis up to the present time has been concerned -with 
the inner or rotating aberrascope structure. 

For analysis, the rotating aberrascope structure is assumed to consist 
of a central cylindrical section and two conical sections. Tho main bearing 
is located midway between the ends of the cylindrical section. The primary 
mirrors are supported at the ends of the cylindrical section. The conical 
sections, which are fastened at their large ends to of the cylindri¬ 


cal section, support the secondary mirrors. This structural arrangement is 
shown schematically in Figure 3^ • 



Figure 3^ 

Aberrascope Structural Arrangement 
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In the structural analysis it is assumed that both the cylindrical 
section and the conical sections will be made of invar. 

a. Structural Deflections Under Static Forces 

Maintenance of the alignment of the secondary mirrors with 
respect to the primary mirrors is the most critical of the structural re¬ 
quirements. lhe aberrascope structure may be regarded as a beam supported at 
the center. Hie action of static forces will cause tius structure to undergo 
deflections with respect to its support. The magnitudes of the deflection and 
slope at the small end of the conical structure relative to the deflection and 
slope at the end of the cylindrical section determine the misalignment of the 
secondary mirror with respect to the primary mirror due to static forces on 
the aberrascope structure. Since the aberrascope structure is symmetrical 
about its bearing, it is necessary to consider only one side of the structure 
in carrying out the analysis or deflections and slopes 

Figure 34 shows the structural dimensions and forces acting on the 
structure which are involved in the analysis of slopes and deflections. 

= total length of the rotating aberrascope structure. 

L = length of the cylindrical section. 

& = length of the conical section. 

R = outside radius of cylindrical section and large end of conical 
section. 

Rj = inside radius of cylindrical section. 

R = mean radius of large end of conical sectio”. 
c 

R„ = outside radius of small end of conical section. 

Cl 

R = mean radius of small end of conical section, 
m 
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thickness of cylindrical section, 
thickness of conical section, 
force due to veight of secondary mirror. 

force due to veight of primary mirror, 
force due to veight of cone. 

Die thickness of the conical section, t, and the thickness of the 
cylindrical section, t , are the parameters -which are varied in the analysis. 

The analysis of deflections and slopes is based on the assumption that 
the aberrascope is in a force field of one g. Diis assumption is made merely 
for convenience in carrying out the analysis. Since the deflections and 
slopes, in the range considered, are linear functions of the loads, the results 
obtained for a one g loading can be used to determine the deflections and 
slopes for the actual force levels which will be encountered by the aberrascope 
once there force levels are knovn. 

In evaluating the results of the analysis in terms of the structural 
accuracy requirements, it will be assumed that the static loads on the 

-P 

aberrascope structure vill be those corresponding to a force field of 10 * g. 

Obtaining deflections and slopes for a 10 c g field from the analytical 

results for a 1 g field involves simply multiplying the analytical results by 
-2 -2 

10 . Die assumption of a 10 g level for the static disturbing forces on 

the aberrascope doling actual operation is believed to be icry conservative. 
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Figure j.' shows £ cantilever beam equivalent of one side of the aberra- 


scope structure of Figure ' -. 



1antilever Beam Equivalent 


Consider first only the cylindrical section, Figure 3 V of length L/2. The 

deflection and slope cf this section at distance l/? from the fixed end, Figure 

36 , is produced by the combination of a force, a moment, and tne cylinder weight. 

The force, F, is a summation of the weight of two mirrors (F , and F „) and 

ml m 2 

the weight of the conical section, The moment effect in turn is a sum of 
the moments due to the mirror, F^, and the weight of tne cone. Evaluation 
of each of these effects yields the following equations, 

Slope: C = FL 2 (l) 

f 8Il 

Deflection; y_ = Fl^ 

F 2511 


lit 


( 2 ) 




Slope: 


ML 

Sec 


( 3 ) 


9 



w 

(5) 


Deflection y = Wl£ (6) 

V ^EI 

Applying the principle of superposition one can now say that the slope and 
deflection for the cylindrical section are 

9 cylinder = Q_ + 9 + 9 (7) 

r ID v 

and 

y cylinder = y p + y m + y y (3) 

Consider now the conical section of length / , figure 34. The conical 
shell may be regarded as a cantilever beam fixed at ’he intersection of the 
cylinder and cone. Uhe slope and deflection for the loaded conical shell 
cantilever beam, Figure 42 are calculated by the double integration method 
which considers the position of the beams neutral elastic curve. Effects 
due to the concentrated force and the beam weight are calculated separately 
and then combined using the principle of superposition. 
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Figure 36 
Conical Shell 


From Figure t 3 6 the following relationships hold: 


K = k 


Z tan 


ft 1 


R = R 
n a 


- t 

2 


R = R 
c o 


t 

p 


tan j3 


R - R 
c n 


L 


Ohe differential equation of the elastic curve is 


d 2 y _ M 

, 2 " El 

dx 


Where: 


M = Bending moment 
E = Modulus of elasticity 
I = Moment of inertia 



The moment of inertia for a thin walled conical section is derived from 
the configuration of Figure 37* 

Y 



Figure 37 

Moment of Inertia Consideration 

Using the classical definition for moment of inertia one has the equation 


I = £ dAy 2 = j y 2 d£ (14). 

Referring to Figure 37, 


dA = 

dst 

('15) 

dfi s; 

Rdo<, 

(16) 

y 2 . 

(R sino{ ) 2 - 

(17) 


Substituting equations (15), (16), and (17) in equation (l f :) yields: 


J I* a ) 2 (R 6in«) 2 Rdo(t (18) 
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Solving equation (IS) gives: 


1 

£ 


*x = 


t R- 



t R J 



Therefore, the moment of inertia I for a conical section is: 

" irtH 3 


(W) 


Referring to Figure 36 and considering the slope and deflection for the 
concentrated load yields the equation y " = ~ and therefore 


where 


or 


y 


11 


-FZ _ 

TT t e ^F n ~+ z 



- ? Z ■ B » bending moment 
I » + Z tan^J 3 -fl-t 

y 1 «* —~.JL ( z 6. z _ 

7Tt E / (,R U + Z tan/3 ) 


( 20 ) 


( 21 ) 


Integration of equation (21) gives the following: 


y’ 


- -F 
'jrt s 


_1_ L _J_ 

tan/} 1 F = + S *** 


+ R n_+C ’ 

2(R n + Z tasiyfj ) 2 . 
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Solving equation (22) for C, where the conditions apply that y = 0 and 
(Sq + Z tan ) = R c when Z = ,/ c , one contains 


- F 


f- -i 

l— R- tai 


*n 


*** E L R c tan"y5 2 R c ta§ 


+ C 


] 


( 23 ) 


or 


C = 


Rn 

R c tan 2 ^ " 2 R c 2 tan*# 

Substituting the value of C in equation (22) yields 

v' = f- 1 _+ R n _1 + 1 _ 

E jj™ 2 /# \ Ra + Z tan^ 2 (R n + z tan^) 2 J R c tan 2 , 


Rn 


tan 2 ^ 2 R c 2 tan^J 
(24) 


Solving equation (24) for the special case, where y' = e CF and (R n + Z tan/#) 
- Rjj when Z = 0, yields 


/ 2 


" F TTt E (R. 


c_ "_^n 

- R q ) 2 2 Rn' 


P “7 
I 


_i + _i . 

B 2 Rn. R c 2 r 2 


c_i 


or after simplification, the slope, at end of beam due to a concentrated 
load at point Z = 0 is 


9, 


_ / 2 
~ F (4±~) 


•F 2 7TE t Rn ' R 
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Integrating equation (24) for beam deflection yields 


- F 

f Rp 

f- 1 1 

1 1 r/n (Rn +Z tan#) ) 

E7Tt i 

1 2 tan 2 /5 ( 

) tan,<3 (Rn + Z tan#) j 

, tan 2 # ( taQ # j 


+-L__ -____ + c, 7 (26) 

R c tan2^ 2 R c 2 tan^J J 

Solving equation (26) for Cp where the conditions si-:'\ chat y = 0 and 
(Rq + Z tan j3 ) = R c when Z => J c , one obtains 

0 * ■ — ^ f* - ^n __ -/a R c + J c - Rn c. + C~f 

E L 2 ta “ 3 /? ( R c) tan 3 # R c tan^# 2 R c 2 tan^<3 J 

(27) 

or 

c s R n + /n_ R c _ Jc + Rq A 

1 ~ 2 R c tan3^ tan^ R c tan^# 2 P. c 2 tan^5 


Substituting the value of Cj_ in Equation (26) yields 

Rn /n ( R n * z 


-F 


E I- 2 tan ^/? (Rn + Z tan j$) tan 3 # 

__+_ h _+ Ajk. . . A + Rp ^ 1 

2 R c 2 tan^3 2 R c tan 3# tan 3 # R c tan 2 # 2 R c 2 tan 2 # J 


R c tan 2 /? 


Solving the deflection equation (28) for the special case, where y s = y c 

F 

and (R n + Z tan jQ) when when Z = 0, yields 


13.9 




y cp = 


-F 


T_ *n 


/n*n 


‘■•a 


/r 

xi n K c 


E7>t ^ 2 tanS^ (rJ tan3^ 2 R c tan^ tan3# 


/. 




~-i—~ + / 

R c tan'^F" 2 R c 2 tun^ J 


or after simplification, tiie deflection at and of team, due to a con¬ 
centrated load at point Z = 0, is • 




= 2 nt 5 (Rc - R n )3 






(29) 


Now considering the slope and deflection due to the team weight. Figure 3^> 
one again has the equation 


■> u 

a y M 

dx 2 = EI 

where E and I are as in equation (20) and M is a function of weight and 
distance. 

The moment equation is 

M=/?^ZAdZ ( 30 ) 


Taking any distance Z one finds that the cross-sectional area. A, is given 
by the expression 


A » 7T(R + §)* - 7T(R - J) 2 
or 

A = 2 77R t 


( 31 ) 
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Evaluating R, in equation (31) j as a function of Z gives 


A = 2Tf'(R + Z tan P ) t 
n 

Substituting equation (32) in equation (30) and integrating gives the 
moment as 


(R + Z tan/3) ^ R (R + Z tan #) 2 

M = 21ft f - 5 —-2—^- 5 -- + C f 

3 tan j8 2 tan ft 

Solving equation (33) for Cq where M = 0 when Z = 0 gives 


u 0 - - x— - -p—r v. 

2 tan ft 3 tan ft 

Substituting for C Q in equation (33) gives the moment equation as 

(r - Z tan #) 3 R (R * Z tan /?) 2 R 3 R 3 

M = 2T\f > t -2- —- - -2- + - 2 “2-(35) 

3 tan ft 2 tan ft 2 tan ft 3 tan ft 

Substituting equations ( 19 ) and ^ 35) in equation ( 13 ) yields the equation 

y „ _ g| ft /° E(”, »^ ) 3 - Vl 

6 tat [R r + Z tar^j 3 E 
Integrating expression ( 36 ) yields the slope equation 


y , = _/L_ 

3 £ tan 2 # 




+ Z tan/9) 


5 + C 2 



121 


Solving the slope equation (37) for Cg, where the conditions apply that 
y* = 0 and (Rj, + Z tanfl ) = Rg when x = J c , one obtains 


3 % log R c 
Cg = ---- + 




tan jQ 2 R c 2 ran /3 


l 


( 38 ) 


Replacing C2 with the solved value gives the equation 


, r , „ ? To 7 3 R n log (Fq + Z ton /3 ) 

J ** « ■ !A - I •* " “ ■ ■ — — . . . . . 1 


O 

T> J 


3 E tan d /$ 


tan /^ 2 bzu/S (r q + Z tan ) 


L 


3 Rn log R c + V - 2 / 1 

tan /f 2 R c 2 tan^ " ° J 


(39) 


Solving equation (39) for the special case where y = 9 C and (Rq + Z tan Q ) * 
when Z n 0 yields the slope equation at the end of the beam due to its own 
weight. 




3 E (B c - R n 


55 [3*tj£>♦ | 


(to) 


Integrating equation (39) for beam deflection yields. 


[f ■ ( los H " + 2 jp» 

- - fp 3 — /*- 1 _| 

2 tan/? ’ tan^/3(Rn + Z tan£ J 


+ Z tan ) - Rq log (R^j + Z tan/? ) 

3 Rn Z log R c 


O 

Rn J 


taa^> 2 R c 2 tan^ 


Jz + 


•] 


(41) 
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Solving the deflection equation (4l) for C3, where the conditions apply 
that y s 0 and (r q + Z tan j& ) = R c when Z = c, one obtains 


C 3 = ’ 


> 

^ ^ 108 “■=' ^ {«'' 108 “=}]j 

Rn 3 _ 3 Rn ^/c lo S R c - R n 3 /c + Jz 

2 R c tan^3 tan^ 2 R 0 2 tan 


(42) 


Substituting this value of C3 in equation (4l) yields the general deflection 


P 


3 E tan^5* 


C* 


, 2 - 3 R a 


tan jft 


j^Z log (Rjj + Z tanyS ) - tan^jl 


tan^ 


{(Rn + Z tan 0 ) - R^ log (IL + Z R ))1 ) + -5--- 

( * * * UJ 2 tan 2 ;? (R n + Z tan# ) 


+ 3 > zl ; 8Rc + -^~ - 2 y,-4 2 + i^(/ c io 6 R C 

tan/j’ 2 R c 2 tan$ zanp r 


- tan^ (S 0 " Kn l0 6 R c) 


^ i/ 


. .&L. . R c 

2 R 0 tan 2 /y tan^J 


«n 3 / c 


2 R c 2 tan^ 


r. + 2 , 4 2 i 


(43) 


Solving the general deflection equation (43) for the apodal case, where 
y ■ y c and (Rq + Z tan^) ■ when Z ■ 0, yields the deflection equation 
for the deflection, at end of beam, due to the beam weight. 
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M 4 

y °v = 3 E (R c - Rq)^ 


"c R n 


l + 3 log & - --fe -fe + 3 l 


2 ' -r7"^ ‘ FIT 


(Rc - R, 


, 2 

H) 


n.n 

if*/ 


The deflection and slope at the small end of the conical section due 
to a concentrated force, F^, may be obtained from equations (25) and (2?) 
of ■‘••he above development. 


0 Force = 
cone 


- F 


2 -nr E t R n R, 


(4f 


(25) 


y Force 
cone 


- P Yc 3 


**n ~ R c o /?C\ j. / R c " Rn 


o / \o I R + 2 lo B (r“) •+ (— s - 

2 IT t E (R c - Rq)3 j_ R c Rn Rc 


)( |- 2 >] 


(29) 


The deflection and slope at the small end of the conical section due to the 
weight of the conical section is given by equations (40) and (44). 


9 Cone _ 

Weight = 3 


/°4 3 r 

s (R c - Rn) 3 L 


3 Rq log — + 


*n ' 2 R c 2 2 


| «n - 2 R. 


•] 


(40) 


y Cone 




Weight 3 (r c _ Rn ) 4 E 


it* 


4 1 Inn ( R C\ 7 R C Rn , , ( R c " "U* 


Equations (25), (29), (40), and (44) have been restated on this page 
in order to summarize the results of the previous analysis. 
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The principle of superposition is now applied to give the slope and deflection, 
for the conical section, as 


6 Cone = 6p Cone + Cone (4?) 

y Cone = y F Cone + y w Cone (46) 


Total deflection of the free end on the aherrascope 



is now 


obtained by translating the cylinder deflection to the free end and adding 


cylinder and cone components directly. 



Figure 38 

Beam Problem Parameters 


Translated cylinder deflection, y , (Figure 38) is given by the equation 


where: 


y c » 7 + 4 e i 

1 o. „ 

* c ± -l 


(47) 
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The total deflection of the aberrascope section, -L, 


can now he stated as 


y t = y a + y c (W) 

where: y a = cone deflection 

y c = cylinder deflection 

and the slope of the configuration can he stated as 

©t - ©i + ©p 

where: ©^ « slope due to cylinder 

©2 = slope due to cone 


Utilization of the above equations will provide information as to the most 
desirable cylinder and cone wall thickness consistent with allowable 
secondary mirror deviation. 

b. Evaluation of Deflections and Slope s 

Calculations of deflections and slopes for SDciTSSCOpG 
structure were made using the following numerical values. 


4- 

= 


*VL 3 


5 inches 

13 inches 

J ! inches 

(4 - t c ) inches 

2.5 inches 

(2.5 - t) inches 


t^ = .1, .2, *3, *4, *5 inches 
t = .05,. .10, .15, .20, .25 inches 

*Bl » *°3 # 

- 1.51 # 
f = .294 #/in. 3 
E = 20 x 10 6 #/in. 2 



The data obtained, assuming a one g loading, is tabulated in Table VIII. 

Figure 39 shows the total deflection at the small end of the aberrascope 

-2 

for a 10 g loading. Since the deflection at the end of the cylindrical 

section is negligible compared to the total deflection, the curves of 

Figure 39 also show the deflection of the secondary mirror relative to the 

primary mirror for various values of wall thickness for tne cylindrical 

and conical sections under the assumed loading. Figure s..ows uhe total 

-2 

slope at tne small end of the aberrascope structure for a 10 g loading. 

Tne slope at the secondary mirror relative to tne position of the primary 
mirror is snown in Figure hi. 

c. Weight Calculations 

In addition to meeting the accuracy requirements of the 
instrument, the aberrascope structure must be kept within reasonable weight 
limits- Calculated weights of the aberrascope structure for various com¬ 
binations of cone and cylinder wall thickness are tabulated in Table VIII 
and are plotted in Figure hi. The calculated values include the weight 
of one cylindrical section and two conical sections. Weights of the bearings 
or components mounted on tr structure are not included. 

d. Conclusions from Static load Analysis and Weighs Calculations 
The curves of Figures 39> 40, and 4l show that both deflec¬ 
tion and slope increase with increasing cone wall thickness. Deflection 
and slope decrease with increasing wall thickness of the cylindrical 
section. Therefore, the best structural arrangement for minim-sing mis¬ 
alignment between the primary and secondary mirrors is one having thin 
walled conical section and a thick walled cylindrical section. However, 
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Deflection - inches x 10“ 



Figure 39 

DEFLECTION Mr END POSITION OF THE ABERRASCOPE 
UNDER 10-2 g LOADING 
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Cone wall Sickness - Inches 


Figure 2*1 

SLOPS AT SECCOLAEY MIRRCR RELATIVE TO 
PRIMARY MIRROR 









Cylinder Thickness 









the overall weignt of the structure increases with increasing wall thickness 
of the cylindrical section, so that a compromise must he made between 
structural and weight requirements. Also, the cylindrical section must 
be sufficiently rigid to provide support for the aberrascope main bearing. 

From the above consideration wall thicknesses of .05 inches for the 
conical section and .2 inches for the cylindrical seo~'~n appear to be 
reasonable choices. The weight of such a structure would be approximately 
15 pounds as shown in Figure 42. Comparison of the deflection and slope 
magnitude for a loading of 10 g on this structure with the accuracy 
requirements may be made from a consideration of the sketch below. 



x coordinate lies along optical axis, 
y is normal to x. 

© measures rotation with respect to a plane perpendicular to x. 

The quantities b x, by, and [ 9 indicate small displacements relative 
to x, y, 9 coordinate system. 
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Displacements y and o 9 correspond to deflection and slope 
respectively. The allowable errors in alignment during one rotation of 
the aberrascope due to disturbing forces are: 

£ y ^ 2 x 10"5 inches 
r , 

<5 © .01 seconds of arc 

The calculated values of deflection and slope fo; Uic structure and 
loading assumed above are shown in Figures 39 ana 4l. Using a cone wall 
thickness of 0.05 inches and a cylinder wall thickness of 0.20 inches, 

the following values are obtained: 

c -8 

5 y = 8-9 x 10 inches 

6 9 = .0016 seconds of arc 

The deflection is less than the allowable error by a factor of 200 
and the slope is less than the allowable angular misalignment by a factor 
of 6. 

e. Structural Response to Dynamic Loads 

While the above analysis of slope and deflections under 
static loading showed that conical structural sections with very thin walls 
would meet the accuracy requirements with larger factors of safety, it is 
also necessary to investigate the behavior of the structure under vibra¬ 
tional disturbances. In order to avoid instrument errors due to vibration 
it is desirable to have the lowest resonant frequency of tho structure 
significantly higher than the frequency of any disturbing vibration. 
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Since the cylindrical section of the structure will he much more 
rigid than the conical sections, the vibration analysis is concentrated 
on the conical sections. The ways iu which the conical sections may 
vibrate may be classified as: 

1. Cantilever beam vibrations. 

2. Thin shell vibrations. 

3* Torsional vibrations. 

Each of these forms of vibration will be considered in turn below, 
(l) Cantilever Beam Vibration 



A 

I 

A 


f jj 2 - (d - P-tfj 

$ P ' ( d ‘ 2t) J 

Cro$3 sectional area of shell 
Moment of in inertia of shell 


The first mode of vibration as a cantilever beam is calculated by 
the Rayleigh method. Elementary beam theory is used in the calculations. 
This introduces some error in the results due to neglecting transverse 
shear and rotory inertia. However, the frequencies obtained should be 
sufficiently accurate for the purpose under consideration. 


135 



The deflection curve under vibration is assumed to be the same shape 
as the deflection curve for a uniform cross-section cantilever beam under 
a uniform static load. The actual shape of the vibration deflection curve 
will be somewhat different than the assumed curve. This, however, will 
introduce only a small error into the calculated frequency. 

Taking 

X = a (6i 2 X 2 - 4/x 3 + x“> (3) 


where x = mode shape 

a = a constant 
the boundary conditions of the problem are: 


(± 1 ) 

k dx ' 


X = 


(X) = 


x = 


0 

0 

0 

0 


d 2 X 

(EI ^r> 


4 - (El ^r) 

dx dx^ 


= 0 
X -{ 

= 0 
x J 


W) 


The resonant frequency is obtained from the equation: 


6J 



2 Eg 

V 


(\ (d 2 X, 

1 I( ^ _) 

1 — 
J A X* dx 


dx 


o 

resonant frequency 
Modulus of Elasticity 
acceleration of gravity 
density per unit volume 


(5) 


136 





length of shell 

distance from fix end of shell 


d = outside diameter of shell at any cross section 
Djl = outside diameter of large end of shr'* 

= outside diameter of small end of shell 


Letting 


" °2 
/ 


= C 


( 7 ) 


Substituting equation (7) in equation (6) 


d = - CxJ 

Substituting equation (8) in equation (2) and expanding the terms gives 


( 8 ) 


i = 1 

o 


[ Alx+A2 * 2+fl3x3,A “] 


(9) 


where 

t = shell thickness 

= (etDjL C - 3D! 2 C - 4t 2 C) 

Ag = (3DiC 2 - 3tC 2 ) 

A 3 = ° 3 

Au = (Dl 3 - 3 tDjf~ + ht 2 D x - 2t 3 ) 
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Substituting equation (8) in equation (i): 


A = IT t [(Di - t) - Cxj = IT t £(D]_ - t) - X J 


(10) 


J 


1 (^ 2 ' 

I (—p~) ax 

dx 


3 e 2 7 T t / 5 f / 

= -140-L Al * 


. /2 


+ 3 A, ^ + 168 Aa 


(ID 


j a X 2 dx = g2 ^ — f 144 + 584 - 728 tj 


( 12 ) 


Setting equation (11) equal to equation (12) and solving for (j gives: 


Id 


j - 28 A^ 8 A2 3 A^ 168 Ai| - 

6-75 f j ZTXZ 7 _ j 
r 144 + 584 Dp - 728 t J 


(13) 


Equation (13) may be solved for u) • 

(u) Numerical Calculations 
D x = 8 
% - 5 
./ - 13 
c , tl = .231 

E = 20 x 10 6 

y> = .294 
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t 


.2 


f - 

a 2 fr 


(b) Results of Calculations 

t = .05 

f n = 1665 cps 

t = .1 

f n = 1*48 ~ps 

t = .2 

f n = 1616 cps 


The resonant frequencies for the conical shells for cantilever 
beam vibrations are verj' high. Variation of the thickness over the 
range considered makes little difference in the resonant frequency 
(about 3$). This is probably less than the error in frequency due to 
the various simplifying assumptions made in the analysis. 

(2) Thin Shell Vibration s 

The thin shell vibration of a cone may be composed 
of circumferential modes, axial modes, or a combination of these modes. 
If a section perpendicular to the axis of the sneil is considered, the 
vibration may consist of the periodic movement of a number of stationary 
waves distributed around the circumference. This form of vibration is 
illustrated in Figure 43. The number of circumferential waves is indi¬ 
cated by "S". 
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m = 1 m = 2 m = 3 


Figure 45 . Axial Waves 



Figure 46. Location of Nodes 
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Figure 45 shows the first three axial vibration nodes of a conical 
shell with one frea and one built in end. Figure 46 shows the location 
of nodes for one case of combined circumferential and axial vibration 
(n = 3, m 4). 

The calculation of resonant frequencies of thin walled conical shells 
is a complex and time consuming process. Since the information of 
principal interest is the lowest resonant frequency oi the shell, the 
analysis is limited tc. the case of circumferential vibration only, (m = 0). 
The frequencies of the axial modes rr the combined axial and circumferential 
mode will all be higher than the lowest frequency of purely circumferential 
vibration. 

A Rayleigh-Ritz procedure is used to determine resonant frequencies 
of circumferential modes of vibration of a conical shell which is built 
in at its large ends and is free at its small end. 

The Rayleigh-Ritz procedure equates maximum potential energy to maximum 
kinetic energy. A displacement is, first, assumed for the mode shape of the 
vibration. The assumed mode shape must satisfy the displacement boundary 
conditions and contain arbitrary constants. The maximum potential and 
kinetic energies of the system are obtained from the assumed mode shapes. 

The resonant frequencies are calculated by applying the Rayleigh-Ritz 
procedure to the resulting statement of Hamilton's principle. 

The analytical development follows that given in the paper "Vibrations 

PI 

of Conical Shells" by H. Saunders and E. J. Wisniewski, except that 
modificacions are introduced due to the different boundary conditions of 
the present problem. 



The middle surface strains of the shell are expressed in terms of 
middle surface displacements up Ug, and ^ as shown in Figure 47. Each 
displacement is assumed to be a sinusoidal function of time. 

UT (Zj 9, t) = ur (Z, 9) sin tut 

u^ (Z, 9, t) = u^Z, 9) sin w t 

u 2 (Z, t) = U2(Z, 9) sin ui t 

I 


Figure 4-7. Middle Surface Displacements 

Z = cylindrical coordinate measured along axis of cone from apex 

u^ a meridional displacement at position Z, 9, and time t 

ug a Tangential displacement at position Z, 9, and time t 

U) = inward displacement normal to the tangent plane to the shell 

surface at position Z, 9, and time t 

jIq = Z coordinate to smaller edge of conical shell 

il = Z coordinate to larger edge of conical shell 

q a one-half the included angle of the cone 

h = shell thickness 
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Linearized expressions given by Love 22 are used for the middle surface 
strains, £ p £g, £12 and the < * an 8 es °f curvature, Kp, Kq, and 



£ 2 = .. U 1 cos<x _ ur cos 2 <x 

2 ^zsin<x + T v '° ^ T “H5o“ 


f _ z cos Oi * (* ) + i c - 030 < -b: 

C 12 - Z cos o' -g— (—) + z sinoC -J©“ 


Ki - ^ C ° S « 


^ u 2 cos ^cx + cos 2 PC ^ 2 uJ . cos 2 o< h ur 
<* 9 Z 2 sin o( Z 2 sin 2 (* c )© 2 z ~ Tz ~ 


,, ^ °2 cos^OC _ cos3 Q( 


C< j 2 gr _ cos 2 ^ o 


^ Z Z sn 0{ g 2 sin OC ^ ^ sinoc ^ Zq 9 Z 2 s -^ noc ^ ® 


Ci = meridional strain of the middle surface of the shell 
£ 2 = circumferential strain of the riddle s’^-fac" 

= shear strain of middle surface 

= change of curvature of middle surface in meridional direction 
due to deformation of the middle surface 

Kg = change of curvature of middle surface in circumferential 
direction 

Kj 2 = twist of middle surface 
The maximum kinetic and potential energies are: 

T = £ //>h (u x 2 + Ug 2 + or 2 ) dA (3) 

A ' 
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(4) 


V-i 


2 ^ _ n( 2) j | 12 + " 2 C 1 ”V) ( K i ^ “ K 12 2 'J 


2^1 


+ (£l+6 2 ) " 2 (?. -7) (£ 2 t 2 - jr i 12 )j dA 


Substituting Equations (1) into Equations (3) and (4) .jgetner with 
proper limits gives: 


1 2 
T = ~ fi h ui 


2rr , .2 2 2 v Z sinOC . _ ... 

■ (uJ + u. + u 0 ) - 7~T~ d Z d© 

. . 1 d cos c 0<. 


(5) 


1 Eh 

V= 2 (1-v 2 ) 


2if 2 g 

. ^ ^(K X + Kg) - 2 (1 - v)(K x Kg - Kj/j 


,7 ^ 


( 6 ) 


+ (t x + e 2 ) 2 - 2 (1 -y)( t 1 t 2 - j^): d Z d 

) COS 2 0t 


T = maximum kinetic energy of the vibrating shell 
V = maximum potential energy of the vibrating shell 
E = Modulus of elasticity 
•y = Poisson's ratio 
& = density of shell material 
uJ = frequency 

The displacement boundary conditions for a conical shell built in on 
the edge at the large end and free at sma)l end are: 


144 



(7) 


{jf = 0 at Z = 


l 


-^LisT = 0 at Z = / 

Tz 1 


The mode shape assumed for or is: 

n2 


or = ( / 1 - zf [_A + B Zj sin S 


9 


(S) 


S = number of complete waves around circumference 
A, B = independent constants 

The other two displacements u^ and ug are selected so that £ and £2 
are zero. 


Then: 


= 0 


Ug = S2|2L_ ( /l - z f [a + BzJ cos 


S 6 


(9) 


'12 = 




cos SG 


Ki = cos 2 o( |^2A + (6 - 4)j sin S9 

[ r /), 2 j/y 2 . . -“ii. / r rt , 1 

|A j^/?(l “ + 2 (l g~; + 


Kg = COS 0( 


4, 


u n \ 

+ - 5 ^) + {~£- - h + —i)| sin S9 

z 71 z J 


,, ^3cos 2 c( sin 2 a. r . ! A 2 , \ ^ „ / /, Z \ 

K 12-if sin^T ■ L ^ ( 77 _ 


cos S9 


Substituting Up ug, and or from equations (9) and ( 8 ) into equation ( 5 ) 
and carrying out the integration gives the maximum kinetic energy. 
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The corresponding expression for the potential energy is obtained by 
substituting Kp Kg, and £ 12 from equations (9) in equation (6) 
and integrating. 

For convenience the potential energy is divided into two parts 

V « v b + V e (11) 


V b = the portion of the potential energy due to K]., Kq, and 
= the portion of the potential energy due to E^g- 


E h_ tan (X. 
b = 24 (1 -y 2 ) 


3 

TT cos J o< 
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( 12 ) 


tO|M 


+ 6 log - k (1 - /°) + l» (i - 4--) -4 (1 - It 


/o' 2 


p] 

* o 


*4 j} (1 - -4 > - 3 f 1 - -If) * l/n -4 


/ J 

T + (1 "7 i) 


* 2 U - 7*) - 8 (1 

■<1 4 1 


4^) + U log 4 . 2 -v f (gF?l - ^ 


— 7 - - «-y ' 7 - 5 -; 

/o ^ 4 


** (l-) + 2 log -4 + 2 (1 - —) - 4 ( 1 - ^ 0 ) j 


2 / 2 

+ 2 (1 -V) sin 2 0( A- [~§ (1 - -d^) - 2 log JLk 

s L A 4 

4l\ J y 7 3 / 2 

(1 ' 7 W + A B A < 8 < 1 - 4%) - 6 (1 - ) 

+ y 9 2 fi (1 - -^3 ) - 4 (1 - ) + 12 (1 - _4 ) . 8 log A 

^ 4 A /l /o 

~ 2{1 'To + ^[3 (1 " “4s > " 16 - “4? ) + 36 (1 - 4° J 

4 /1 A 


- 16 log 


A ,1 


7 '2 <1 - ~± )| * t (1 - _£j ) - 1 M 1 - ) 

/o 40 J /1 /l 


+ 20 (1 - -£2 


A> ■ kloe( T ) * 2V H 


/3 

3 ( 1?’ 


1^7 





/ 2 / 

- 10 (1 - J° ) + 16 (1 - Jo ) . 4 log _A 1 + 6 (i - _4; 

A A y 0 J 4 5 


- 14 (1 - 


/a / i 

-r-r ) + 10 (1 - -A 2 > + 2 (1 ) sio 2 <y n a . 

4 1 4 J S 2 [3 


A 3 . , 4 2 y , _ x 

—j) - ift.-Sj, - Ml .^2 )ti: , c . _£i \ 

/1 A A /, / 


*P 2j2 i( 9(1 - /* - 16 (1 '"Tf> - !(>-4) 

' 'i A A 

+ A Jjj (l - —2 jj- ) - ^ (1--2^ ) + 3 (l - 4^- ) - 4 (l - ,_4 ) 


/ J 4 '3 • 2 

log --p + 2 y3|| (i _ ) - ~ (1 - ) + 6 (1 - ) 


r, A 


- 6 (1 - *y ) + log -A 

A An 


-rl + ? -mi- A 


/ 3 

' n 


j 2 

+ 11 (1 - y“2 ) " 8 (1 “ "T ) + log -j- + 2 vjjQ j| 


(1 - A_ ) 


t (1 “~/3 ) + 7 ^--4?) -4d- A } 7 + |(i- A) 

^ ^ 1 /1 J y/l 


1 






Since the assumed functions u , Ug, and ur are linear sums in the 
independent constants A and B, the expressions for the maximum kinetic 
and potential energies become quadratic forms in A and E» The Rayleigh- 
Rita. procedure applied to Hamilton's principle gives 
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4t (t - v) = o 

d « 


_jL (T - V) = o 
i)B 

Introducing the expressions for T and V from equations (10), (32), and 
( 3 . 3 ) into equations ( 1 ^) and carrying out the aixfei- -nation gives: 


^ - -I- /•»-* A 2 (5^-mi* 

a a 


2 3 b 

cos (X \ ( OA 1 f i • o 

■y-’ r z ~7i 


- 1 a - JL ) * | ( 1 - 


, 4 / o 

/ o \ , it /-> _ -/ °- ) 

77 3 A 3 

J'l 

1 - / o \ 2 /, 


2 ^ ;2 ) + 2 7 ^ "TT ^ " | t 1 " 1 ~7T ) 

/ 1 J /l /l 


+ | (l+ ° )-(l 

/ O 


4\ 

77' 


1 3 

> * § ( 1 '-73> 

41 ■ 


S h tan tX. 
2b (1 --y 2 ) 


7Tcos 3 oc / l <<^ 2 A <2 (1- -j-£ ) 


2 y 2 / / 

>>0 fi d--A> + ^log + 

L /l / 0 /l 


X, /i £ 


)J + V? 5 (i - Jp ) - 3 (l * ) + 3 IX>S 


(1 - ~ T ~ )1 + 2 (1 “ -72 ) - 8 (1 - - j - ) + 1* log -A 

< ° J / x / 1 / o 
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+ 2 y 


fM 


/2 / / / 2 
°) -»a - -4?-) + 2 logj ♦ (i. z°.) 


A 


-Mi — 

A 1 


A J ' A 
2 


-1 a ^ 


/, 


+ 2 (1 - 7 ) sin 2 c( -A.. f | (1 - j^g) 2 log — 

s *- J 1 •' 

\ / 3 / ^ 

>]) * <Mi 11 --ys) 


. 2 r o 


4 - / *■ / / 

V* " j | U - —A) - Mi - ~^) + 12 (l - -1£) - 8 log _Li 
L A Yi /i / o 


- 2 


/ _ / 3 / 2 a 

(1 - J +/?[t <! - — 3 ) - 16 (l - + 36 (1 - 


/ / -i / 3 / 2 

16 log — 1 - 2 (1 - -4 1 .) + 4 (1 - _J_£_ ) - 14 (1 - _1£_ ) 

Yo J / 3 / 2 

*• X * 1 


/< O 


20 (1 - -^2. ) - 4 log - 2 -v ® (1 - -4L° ) - 10 (1 ) 

Y1 ' ° ( /i /] ■ 


/. 


/ 2 

4 o v 


A, ... A 


/ 3 / 2 

+ 16 (l - JL£) - 4 log _ii 1 + 6 (l - - 12 _) - 14 (i L 2 _ ) 

/1 /o J 


A- 


A 2 


+ 10(l-_£)j „ 2 (1 . y ) .in 2 * (1-4^) 
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;2 ^ / 3 / 2 

d-- 4 %) + ^ 2 |f 

/Ti L. A /1 


+ 12 (1 - J° ) - 8 log h - 2 (1 - A )~j +yS jf (1 - 4% ) 

/l /c /° J L /l 


- 16 (1 - , A-- ) + 36 (1 - ) - 16 log —) - 2 (1 - -A )1 

J 2 71 /o /oj 


/ 3 

+ 4 (1 - JjL ) - .4 (1 

4 3 


+ 2~y 


M 


/ 3 

8 (1 - ) 
3 


4% > + 20 (1 - J° ) - 4 log A 
Al A * 0 


10 (1 - -4^ ) + 16 (1 - -IS- ) 

A A 


J /3 /2 / 

- 4 log 1 + 6 (1 - -12- ) - 14 (I - -4° ) + 10 (1 - -±£ 

AJ 4 3 A 8 /i 


+ 2 (1 - -y ) sin 2 0( 


r / 3 d 

[l 


(1 ' "A" ) + 4 log - 7 ^]) + 28 A 2 ( 9 (1 " -j\ } 


( 1 --^) + 8 (l--^ ) +/ 3 2 jj ( l.-i 2 J . 
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/ 3 / 2 j 

~ 3 ^ 1 73 ^ + 3 t 1 -^ ) ” 4 (l - —) + log 

■h /l /1 


) + log M 


r y 4 / 3 

; 3 (1 j 3 


7 2 3 ) 

A /1 


- 6 (1-r-2 ) + log __4 j + 9 (! . _j/o_ ) _ 8 (l - _A_ ) 

^ A 1 ^ 1 


,n|i -| M ( i -| ) +io8 7; +2 "p[i (i -^> 


3 ^ " y3 ^ + 7 (1 - -~A ) - 4 (l - -^4 )|+ | (1 - _4_ ) 


/ 3 / 2 p 

12 (1 - -±^ ) + 11 (1 - Jo j ” 4 (l - —L£ ) 
/l /i /4 


+ 2 (1 - 7 ) 8ln 2 0( j^2 (l - -~ 


4% s. 4 3 . 

1 1 A 1 


+ (1 - 


0 1 \\\ E h TT tan (X cob^C^^i / / ’ n J ' 

T? >] ))) ■ ~ 7) S 2 f A A[i (X - 42, 
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) 


1 

3 


(1 - 





) 



r 





/< 


A 


+ l^ 



( 16 ) 


The two linear homogeneous algebraic equations (15) and. (16) may be 

2 

solved simultaneously for the square of the resonant frequency, (jj • 

J o = 21.666 

J ± = 34.666 

0 £ =■ tan -1 (.115) 

^0 = 7.61 x 10" 4 
E = 20 x 10 6 

-y = .29 

Substituting the above numerical values into equations (15) and (16) gives: 


= [A] [MSB * 10 2 (h 0> 2 ) (^J^) 

- (h 3 ) 8.93044 X 10 8 - 9.97709 x 10 6 (■ :2 75i^—) 

2 c. -986 - S 2 2 

+ 3.62061 x 10 6 C? 88 ' —) + 2.06560 x 10 6 ( _;£LJ-) 
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- 2.87160 x 10 11 ^ bJ J"!!.57616 x 10 3 (h w 2 ) 


- (h 3 ) <^1.53066 x 10 10 - 3.49933 x 10 8 (^936 _-_s£\ 


2 7 -986 - S 2 p \ 

+ 8.60880 x 10 7 (:&°6 V s .. ) + 4.18490 x 10 ( -013 ^ j) 

S~ w 


- 5.98823 x 10 12 ( 


?>] * 0 


= JTaJ £ 11.57618 x 10 3 (h W 2 ) ( s2 ) 

- (h 3 ) ^ 1.53067 x 10 10 - 3.49983 x 10 8 (-- 98 ^ s2 ) 

7 22 .986 - s 2 . 

- 8.60880 x 10 7 (■ ^ - 88 1 - : s ) + 4.16490 10 7 ( * 01 3 _ ) 2 \ 

• U18 s 2 / 

- 5.98823 X 10 12 (-£-)J + ^ B ] [2.96673 X 105 (h uu ? ) 

- ( 4 3 ) / 3.80788 x 10 11 - 1.14657 x 10 10 ) 

+ 1.9993 X 109 — ) 2 + 8.64814 ^ S 2 ) 2 ^ 


1.28199 x 10 lU (_^_) z= 0 

s 
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Substituting values of S and h into equations ( 17 ) and (18) and setting 
the determinant of the equation equal to zero will yield the equations 
which may be solved for the resonant frequencies,(</ . Thus, in order to 
obtain a preliminary estimate of the resonant frequencies of a conical 
shell, an analysis based on the vibration of a thin walled cylindrical 
shell was made 

The resonant frequencies of a conical shell having a small semi¬ 
vertical angle and a small length-to-radius ratio difler only slightly 
from the resonant frequencies of a cylindrical shell whose diameter is 
equal to the mean diameter of the conical shell. A freely supported 
cylinder is considered where the ends are forced to remain circular. 

The slope in the axial direction is not constrained to be zero as it is 
in the case of a fixed end. 

The end condition at the large end of the aberrascope conical shell 
structure will be somewhere between a fixed end and freely supported end 
condition, since the fastening between the aberrascope cylindrical structure 
and the conical structure will not be completely rigid. The end condition 
at the small end of the conical structure will depend upon whether or not 
a stiffening flange 01 ring is used ai. the end of the shell. If a stiffening 
ring is used, the end condition will be between that of a fixed end and 
a freely supported end. If a stiffening flange or ring is not used, the 
small end of the conical structure will be a free end. Thus, the end 
conditions taken for the mean diameter cylindrical shell will not corre¬ 
spond exactly to the end conditions of the conical shell in either case. 
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In the case of a stiffening flange or ring at the small end of the conical 
shell, the assumed end conditions are conservative. If the small end of 
the conical shell is not stiffened, the actual vibration frequencies will 
probably be lower than those calculated for a mean diameter cylinder. 
However, the frequencies obtained from consideration of a mean diameter 
cylindrical shell with freely supported ends should give satisfactory 
approximations to actual conical shell frequencies to serve as a guide 
for the selection of conical shell design parameters for more exact 
analysis. 

The equations used in the calculation of resonant frequencies of the 
mean diameter cylindrical shell are: 


f - 2 p a V~a~~ 1 

f a frequency - cycles per Becond 
a » mean radius of cylinder 
y° « density of cylinder material 
E ■ modulus of elasticity 
■y * Poisson’s ratio 

A *» frequency factor obtained from Figures 4, 5, 6, and 7 
of reference paper. 


A 


me an circumference 
“ uiial wavelength 


mlTa 



m = number of axial half-waves 
y » length of cylinder 


thickness t 
maan radius “ a 


t b thickness of cylinder wall 


(?) 
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The calculation of resonant frequencies was effected for a cylinder 
fabricated of Invar. The significant parameters of the material and 
dimensions of the cylinder are listed below. 

E = 20 x 10 6 
yO = .294 

a = 3.25 

g = 386.4 

t = . 05 , . 10 , J. 5 , . 20 , .25 

J = 13 


Table IX 


t 

.05 

.1 

.15 

.20 

.25 


.015 

.031 

.046 

CM 

VD 

O 

.077 


Table X 


m 

1 

2 

3 

4 

5 

m 

,785 

1.57 

2.36 

3.1^ 

3-93 


Calculations were carried out for n = 2 , 3, 4 , 5 and m = 0 . 1 . Larger 
values of n and m would give frequencies which are too high to require 
consideration. 

The resu’ts of the calculations are shown in Table XI and Figures 4 G 
and 49 . Only those frequencies at integral values of n la Figures 48 and 
49 have any significance. The dotted curves merely interconnect points 
corresponding to the same shell thickness. From Table XI, the lowest 
naturax frequency is seen to be 108 cps, which is comparable to the motor 
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Frequency - Cycles per Second 



Number of Circumferential Waves - n 


Figure 48 


FREQUENCY OF VIBRATION WITH NO AXIAL WAVES 
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F requerr-.y 


cf Circumferential Waves - n 
Figure 49 

.f Vibration vith One Axial Wave 
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Table XI 

























the other was made. A value of 1,985 cpa was obtained. This indicates 
that the lowest resonant frequency of torsional vibration is high enough 
so that it will not be excited in the satellite environment. 

f. Thermal Effects 

Some discussion of structural deflections due to thermal 
effects was presented in the quarterly report of September 19. I960. Ho 
further work has been done on this problem. 

g. gyroscopic Torque 

IDuring satellite maneuvers, gyroscopic torques can be 
generated by a rotating aberrascope structure. These torques could cause 
deflections of the aberrascope structure or difficulty in satellite at¬ 
titude control. The magnitude of the gyroscopic torque is given by the 
equation: 

T = I A W A UJ S 

where: 


I a moment of inertia of rotating rVerxascope about its 
axis of rotation 

OJ = aberrascope angular velocity 

(p a component of satellite angular velocity about an axis 

8 perpendicular to the aberrascope axis of rotation 




f “cyl. 


r 2 

L*o + 


. + . 2 -| ( ^2 (y* 3mi [(Rp-t)* - (R s -t)^J 

(Rc-) J + [ 10 (^3 . ^5) 10 [(Ro-t) 3 - ( R a ~ t ) 3 J 


Values of I„ corresponding tc several different combinations of 
cylinder and cone wall thickness are shown in Figure 50* Using these 
values of I a and U) a equal 90 revolutions per minute, gyroscopic torques 
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for values ofU7 from 10 to 60 minutes of arc per second where calculated, 
s 

The torque magnitudes obtained are plotted in Figure 51* To minimize 
gyroscopic torque effects, structure wall thicknesses should be kept as 
thin as possible consistent with structural requirements, 
h. Summary and Conclusions 

For purposes of analysis, a structure consisting of a 
cylindrical section and two conical sections was assumed for the aberra- 
scope. It is believed that this assumed structure is a good model of 
the most practical design for an actual aberrascope structure. Analysis 
of the response of the structure to static and dynamic loading for 
several different wall thicknesses of both the cylindrical and conical 
sections was made. 

From consideration of static deflections, cantilever beam vibration, 
and torsion vibration, the wall thickness variations of the conical sections 
over the range of values considered were not significant. The analysis 
indicated that conical sections with a .050 inch wall thickness would 
meet the structural requirements relative to suable loads with substantial 
factors of safety. A wall thickness of .2 inches for the cylindrical 
section should meet the static load structural requirements. From the 
point cf view of minimizing weight and gyroscope torques the walls of 
both the cylindrical and conical sections should be made as thin as Is 
compatible with other structural requirements. Practical limitations on 
wall thickness reduction are set by consideration of thin shell vibration 
support of concentrated loads and mounting attachments, and in the c se of 
the cylindrical section, support of the main bearing. The thickness of 
.05 inches end .2 inches ior the conical and cylindrical sections would 
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give a total structural weight of about 15 pounds. These wall thicknesses 
should also give sufficient structural strength to meet requirements of the 
main bearing support, and support of concentrated loads end mounting attach¬ 
ments. Although the thin shell vibration frequencies for tne conical sections 
have not been evaluated, approximate calculations based on cylindrical tubes 
whose radius was equal to the mean radius of the conical sections gave a 
lowest resonant frequency of about 110 cps for a wall sickness of .05 inches. 
This frequency is probably higher than the frequency of any disturbing force 
to which the aberrascope would be subjected in actual operation. 

In general, the analysis indicates that an aberrascope structure of 
reasonable weight will have sufficient rigidity to maintain required 
instrument alignment requirements when subjected to the mechanical dis¬ 
turbing forces which are likely to be encountered in operation. A more 
detailed investigation of thermal effects should be made, since these may 
present more serious problems than the mechanical disturbing forces, par¬ 
ticularly for an aberrascope whose structure does not rotate. The strength 
and rigidity of the stiucture analyzed will be substantially improved with 
only a slight increase in weight by use of stiffening ribs and flanges. 

In actual design, flanges are necessary for interconnecting the structural 
components. Thus, we are assured of a structure which will be completely 
adequete for the requirements of either a rotating or non-rotating aberru- 
scope. 
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2 . Detail Design. 

On the basis ef the structural analysis given above, a 
preliminary design of the optical system has been made and is shown in 
Figure 52. The mechanical structure is shown to consist of three cast¬ 
ings, two of which will house mirrors and the third is the central housing. 
The primary cell will mount the primary mirror, the principle light baffle 
and will also provide a means of attaching the optical system to the balance 
of the aberraocope. The primary mirror must of course be polished on Its 
outside diameter to minimize the possibility of cracking during assembly. 

The primary mirror will be seated against a flange in the primary cell; 
hcvcver, the mirror will contact the flange only at three points. At 
these points, a film cf aluminum between .003 " .005 inches thick will 
separate the mirror from the flange. On the back side of the mirror, 
the restraint will be ai three points so that 'the force applied to the 
back of the mirror is directly opposite to the support points on the 
front, fy this means, no bending torques will be applied to the mirror. 

The secondary mirror will be mounted in its housing by the same 
principle as the primary mirror. A spacer is shewn cetweec the secondary 
mirror face and the housing flange. 2ba purpose cf this spacer is to 
allow for adjustment cf position of the secondary with respect to the 
primary* The final stipulation on the position of the focal plana is 
that it be located appropriately with respect to the optical system 
mounting flange. In order to hold the dimension from the focal plane 
to the mounting flange to a tolerance of hkD. 005 inches, the position of 
the secondary mirror must be adjusted. This adjustment can be effected 
by changing the length cf the spacer in fr.nt of the secondary and by 
adding (or subtracting) an appropriate spacer on the back side of the 
secondary mirror. 
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The iexicr. .t t'.e 2 e:,>oisry --1L '.?■ trat tie rue. 

requirements Vetter, th* se'-.'isc:- airr t m,..:.+-irg 'izjaetsr sad the cell 

nuuriirg i « :e cat ':e set. 3,. tr- pexrisry to the osmtral 

hewing on $ f^ab face, a minimum rf Ustcrtioc. o?.~\xn when the ^khtenere 

t*. 

are secured in place. 

The central h:using is simply a conical shell which constitutes the 

major portion of the cpticai system housing. The material which will hi 

used for this housing as well as t"e mirrrr cell housju. v > is Invar. For 

the first pr:t"type acdel_, 3 M- re' zmmends the use cf commercial Invar and 

Vyccr brand glass #7300 (Corning Oiass Works). These materials have.. 

closely matched coefficients cf expansion, however the coefficients are 

larger'than the coefficient if expansile of fused silica. In the 
2(5 * 

quarterly report. OMI recommended the use. cf fused silica for the 

mirrors and a 1 -v expansion alley (lA -665 produced by the Osneral Com- 

municaticna Company) for the structure. This low expansion alloy is presently 

being experimented to demonstrate its cast-ability. Should casting prove 

feasible, the combination of IA*665 and fused silica will be used in future 

instruments. Fc-r the present, cccnarcial Invar and tr ycor glass will meet 

the instrument retirements, sin« the extreme erf'drcnmental problems win 

not be encountered in the first prot:type. 

Gffi is not attempting to perform any lens- 'esign oc. the optical system; 

• however, the basic parameters nave esen worked cut and this informs tier will 

in turn be supplied to the people from when GMI will purchase the optics. 

* 

The mechanical components of the optical system will be fabricated "in 
house" and furnished to the optical supplier who must finally mount the r>ver-al 
optical system specification. 

The type of optical system wj.11 be a varimt from the true Cassegrain 
in that a spherical secondary will be used. This- is known as a Dahl-Kirkham 
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system. The effective focal length will be 100 inches + 4 inches with 
an aperture of k inches. This is essentially an f /25 optical system. 

The system tolerances which have been established are the following: 
bulCi'.MllCC uii Wit din tuiitt from the focal plane to the mounting face 
of the optical system will be + 0.005 inches; the diameter of the blur 
circle at the focal plane containing 6 o£ of the energy will be 0.0025 
inches; and the location of the focal point from the optic."! axis in 
the focal plane will be + 0.010 inches. The requirement for a plane 
mirror surface on the back side of the secondary mirror is for the 
purpose of alignment in test of the complete aberrascope: to align the 
optical axis of the aberrascope with the collimators which will simulate 
stars. 

The optical system design is Buch that it can be completely dis¬ 
associated from the balance of the aberrascope; i.e., no electrical 
connections are required in this unit. The two optical systems required 
for one aberrascope will be identical. The problem of designing the 
aberrascope to track a pair of stars which are not exactly diametrically 
opposed can be solved by designing a pair of wedge-shaped spacers. These 
spacers will be inserted between the optical system mounting flange and 
the central portion of the aberrascope which houses the error detector 
and the electronics. With the use of these wedge-shaped spacers, the 
optical systems can be interchangeable and also, only the spacers need be 
changed if the star selection need be changed before the design is 
completed. The machining requirements with this type of design are thus 
relaxed and the tilt between optical systems can be introduced with 
spacers, not with complicated machining of the basic structure. 
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C. Rotating Error Detector 

Die rotating error detector measures the pointing errors of the two 
Cassagraiman optical systems, ihe electronics which is most immediately 
associated with the error detector yields the pointing error in polar 

MJ&* 

coordinates, and additional electronics converts the polar error signal! to 
cai'-tesian coordinates. Ihe cartesian coordinate error r'^nals are used to 
position the respective optical systems, tending to reduce the pointing errors 
to zero. 

1. Principles of Operation 

Light from the Cassagrainian optical systems falls upon the 
rotating knife edge prisms which are driven by a hysteresis motor. See 
Figure 53. 'Die reflected beams from the two surfaces of each prism fall 
upon ellipsoidal mirrors, and are reflected from the surfaces of these 
mirrors into phototubes. Die surfaces of the inner ellipsoidal mirrors 
(those closest to the prisms) are silvered in seivated bands, so that the 
light beams are optically switched alternately into one phototube and then 
the other as the prisms are rotated, 'ihe electrical signals developed by the 
phototubes are then combined with electronic logic to produce the error signals. 

2. Optical System 

'ihe inner ellipsoidal mirrors are silvered ,in serrated' bands as 
shown in Figure 54. Light reflected onto these mirrors by the rotating knife 
edge prism will pass through the mirrors at the clear areas and will be 
reflected- at the Silvered areas. Die two inner mirrors are identical in 
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|_j A Phototube 




B Fhototube 


Figure S3 

SIMPLIFIED DIAGRAM OF ERROR DETECTOR 
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shape and are silvered identically with the exception that the serrated band 
on one mirror is rotated 10 degrees in relation to the other. The relationship 
is shown in Figure 54 . Ibis figure erroneously shows the silvered bands of the 
two mirrors as being concentric in order to illustrate the phase relationship 
between them. 

Light which is reflected from the silvered surfac-.ee of the inner mirrors 
is reflected directly into one phototube. Light which passes through the clear 
areas of the inner mirrors falls upon the wholly silvered outer mirrors, and 
is reflected into the other phototube by another optical path. 

As the prism rotates, four sequential optical situations exist. These 
are illustrated in Figure 55* Bit situation changes for each 10 degrees of 
rotation of the prism, so the sequence is repeated nine times for each 
revolution of the prism. The entire sequence for one revolution is given in 
Table XII. 

3 . Mechanical Configuration 

Hie error detector is contained in a cylindrical housing with 
square end caps. Windows are provided at the center cf each end cap of the 
housing to permit entry of light from the Cassagrainian optical systems into 
the detector housing. Windows along the sides of the detector housing.pass 
light from the ellipsoidal mirrors to the phototubes. Each of the windows 
must be high quality plane-parallel with minimum wedge angle. The field of 
the hysteresis synchronous motor fits inside the cylindrical housing, and the 
rotor Is suspended at each end with an air bearing. Hie reason for selecting 
air bearings is to avoid the possibility of oil fogging the optical surfaces, 
and to avoid mechanical noise. Figure 55c shows the optical lay out in 
greater detail. 
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OPTICAL CONDITIONS IN ERROR QETirCTOR 
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OPTICAL CONDITIONS IN ERROR DETECTOR 






















The air bearings will probably have the shape of a hemisphere to provide the 
rotor with both lateral and thrust support. The shafts which contain the 
knife edge prisms will protrude through the centers of the air bearings. 

The electronic logic circuits which convert the phototube signal a to 
coordinate error signals must be synchronized with the rotation of the 
prisms. To accomplish this, a coating or strip of magnetic material is 
attached to the circumference of the motor rotor, and a ma^etic pickup 
is placed close to the surface. Synchronizing pulses and reset tones are 
recorded on the magnetic material, and this information is repeated in the 
magnetic pickup as the motor rotor rotates. Magnetic shielding is used to 
prevent the magnetic field of the motor from disturbing the synchronizing 
tape. 


TABLE XIX 


Pointing r< Direction 
Side A 


3 

k 

5 

6 

7 

8 


Phototube A Sees 
a + d 
a + c 
b + c 
b + d 
a + d 
a + c 
b + c 
b + d 


33 

3 ^ 

35 

36 


a + d 
a + c 
b + c 
b + d 


Phototube B Sees 
b + c 
b + d 
a + d 
a + c % 
b + c 
b + d 
a + d 
a + c 

• 

b + c 
b + d 
a + d 
a + c 
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Error Detector Logic 

The required polrr error signal at end "a b" is (a-b). Similarly, 


the required polar error signal at end "c d" is (c-d). These signals can be 
obtained by switching in the sequence shown in Table XIII and then smoothing 
to obtain the polar error signal as a function of time. 


TABLE XIII 



Rotor Position 

Error at a,b 

Error at c,d 

- 

1 

A - B 

B - A 

- 

2 

A - B 

A-B 

’ 

3 

B - A 

A-B 


4 

B - A 

B - A 

1 

5 

A-B 

B - A 

1 

6 

A-B 

A-B 

! 

I 

7 

B - A 

A-B 


8 

B - A 

B - A 

< 

!* 

33 

A-B 

B - A 

i 

3^ 

A-B 

A-B 

■ 

35 

B - A 

A-B 


36 

B - A 

B - A 


It is significant that this switching sequence gives the correct error 
signals when bv/sn phototubes are working or if only one phototube is working. 
When both phototubes are working the gain of the error signal iscK» If one 
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phototube becomes inoperative, tbe error signal is reduced to cc / 2 . It 
will be necessary to include the capability of detecting one phototube 
operation and accordingly increase the electronic gain by a factor of 
2 during periods of single tube operation. 

In addition to the desired error signal, there is an undesirable AC 
signal at nine times the frequency of the rotor speed, i.e. polar error 
signal. This undesirable signal can be raised Lu sr.y odd harmonic of 
the rotor speed by changing the number of segments on the inner mirrors, 
but the width of the light beams compared to the segment width limits 
the number of segments which can be used, 
a. Polar Coordinate Logic 

A series of pulses and a burst of tone are recorded on 
the magnetic material attached to the motor rotor. The pulses are 
phased to the transition of mirror segments, and the burst of tone is 
located between pulses 36 and 1 . Die tone ( 20 -hO KC) is an activation 
signal for resetting the detector logic, and the actual reset is accomplished 
by pulse 1 . The error detector logic is reset with each revolution of the 
rotor to prevent loss of synchronization which could he caused by a 
stray or missed synchronizing pulse. Any error from either of these 
causes can exist only for one rotor revolution. 

Pulses from the magnetic picJsup are applied to a 6 stage counter. 

This counter has a capacity of 6 k bits, but only 36 are used. The first 
two stages of the counter [h bits) are used to operate the polar coordinate 
logic. 
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Binary to decimal equivalents of the counter are shown below: 



Binary Code 

00 

01 

10 

11 

100 



Arabic Code 

1 

2 

3 

4 

5 


As shown in the schematic diagram of Figure % a diocfe matrix obtains the 
error signals for end "a, b". When the Flip Flop §2 is in the zero condition 
the erior signal is obtained by taking the output of phototube A minus photo¬ 
tube B. VZhen Flip Flop §2 is in the one condition the error signal is obtained 
by taking the output of phototube B minus phototube A. Stated in other terms, 
binary counts 00 and 01 yield A-B and binary counts x0 and il yield 3-A. Circuits 
used to obtain the polar error signal for end "c, d" are slitjhtiy more compli¬ 
cated but still straight-forward. For end "c, d", binary counts 00 and 11 
yield B-A and binary counts 01 and 10 yield A-B. End "c, d" requires 4 times 
the diode matrix circuits as end "a,b". 

5he filtered output of the error signal as a function of error detector 
position has the general form shown in Figure 57a. ihis wave form can be 
separated into its X and Y components shown in Figure 57 b and c. 
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b. Cartesian Coordinate Logic 

Die error signals b and c are in the form of a suppressed carrier. 

To transform these error signals into cartesian coordinates it is necessary 
to use a synchronous demodulator. Die circuit is simple, requiring only 
two diode bridges and two transformer's. Die reference signals for the demodu¬ 
lator originate in the 36 bit counter. 

Diode matricies are used to determine the follo-.dng counts: 1, 10, 19, 

28. Diagrams of these matrices are shown in Figure 58 through 6l. For 
the "Y" component the diode bridges are switched at counts 1 and 19. For 
the component the diode bridges are switched at 10 and 28. Die same logic 
applies for the error signals at both ends. 

Die processes required to generate the polar and cartesian coordinate error 
signals resulted in significant harmonic generation at 2, 4, 0 , and 27 times 
the rotor rotational frequency. Diese harmonics require filtering, and the 
requirements and effects of this filtering are discussed in the paragraph 
devoted to servo problems. 

Die binary counter will be reset after each revolution of the motor roior. 
Die tone burst will be separated from the synchronizing pulses with an audio 
fixter, and will be rectified and filtered to provide a reset pulse. 

c. Synchronous Demodulation 

A suggested circuit for synchronus demodulation is shown in Figure 62 . 

t 

Die demodulator is a conventional double ring diode bridge which requires two 
transformers per modulator. A transformer with a center tajpai secondary is 
used as the input for the polar error signal. In accordance with standard 
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practice the diodes should be operated far into their conducting region during 
the on condition, perhaps half of maximum current. 

References for the demodulation bridge are taken from a flip flop. The 
signal from the flip flop is transformer coupled from one of the emitter legs. 
It will be noted that the flip flop has two reset inputs, labeled o j6 and 
and no counting input. For the "X" axis the o p input i.s triggered by the-1 
count and the Tifi input is triggered by the IQ count. For the "Y" axis- the 
o j> input is triggered by the 10 count and the fr j> input is triggered by the 
28 count. 

A total of 4 synchronous demodulators will be required for each error 
detector. 

The circuit shown was selected because it illustrates the important 
functions while remaining free of electronic sophistication. It is 
recognized that the final circuit will have many refinements. 
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D Art j cula~ion of Rotating Error Detector 

l- Detector Mounting 

The rotating error detector must be mounted at the center of 
•he aberraseope such that the axis of the rotating knife edge is parallel 
with the optical axis cf the inst-ument, Of the six possible degrees 
of freedom, this rotating error detector must be capable of having 
four, two in translation and two in rotation. Restraint m 's~ be pro¬ 
vided for the other two possible degrees of freedom: rotation about the 
optical axis and translation along the optical axis. The straightforward 
manner in which this can be accomplished is by mounting each end of the 
rotating error detector in a pair of linear slides as in machine tools 
where the vertical slide is carried on the horizontal slides. This, 
of course, is a complicated method of accomplishing the end result and 
involves many sliding surfaces. In a hard vacuum environment, it is 
necessary to keep the number of moving parts to a minimum; thus the 
approach of one slide upon another will be ruled out. 

Consideration was given to the use of linkages to perform the 
required task, but the difficulty here is the cross coupling between the 
horizontal and vertical motions. To separate the two, one mechanism would 
be required tc carry the second mechanism, again involving a large 
number cf parts. The most desirable mechanism is one which could provide 
independent motion. 

Since the final motion required is simply a translation (the rotation 
of the error detector package in a given plane can be effected by difference 
in translation of the two ends), consideration was given to using simply 
four push-reds to articulate the ends of the error detector. Restraint 
from rotat.on about the optical axis could be provided by shaping the ends 






of the push-rod like a "crow's foot". An end view of the error detector with 



The ends of the rotating error detector will be shaped as a square and 
thus the horizontal push-rods shaped as shown in Figure 63 will prevent 
rotation of the entire package. 

The requirement that no rotation or +ho o n tj re package take place 
should not be confused with the requirement that the knife edge rotate 
inside of the package. 

Figure 64 illustrates the side view of the detector mounting. 

Restraint along the optical axis is provided by making the end of one of 
the push-rods a spherical end. This spherical end is seated in a precision 
"V" groove normal to the optical axis. Only one such axial restraint can 
be provided; it is obvious that two such restraints would yield a redundant 
support. Thus, if the vertical drive articulates the error detector 
package, it will simple slide across the faces of the horizontal drive- 
rods. 
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Lubrication of these surfaces can be accomplished by the use of 
molybdenum disulfide. This lubricant is suitable for surfaces where 
sliding contact is taking place. 

2. Detector Drive 

As mentioned above, the error detector package will be 
articulated by four push-rods; two in the vertical plane and two in the 
horizontal plane. The package will be retained in place with four 
additional push-rods which will spring load the detector package against 
the driving push-rods. 

Figure 64 illustrates a drive mechanism which uses servo motors. 

Size 8 components are envisioned. Each motor (four are necessary) will 
have a gearhead and a tachometer, the tachometer will be used for rate 
feedback into the servo loop. The motor will drive one stage of external 
gearing which in turn will rotate a cam. A cam follower roller will be 
attached to the push-rod; thus the rotation of the cam will be transmitted 
into linear motion of the push-rod. The rise of the follower roller will 
be linear with respect to 'the rotation of the cam. The cam center will be 
offset from the push-rod axis in order to provide an advantageous 
pressure angle to the follower roller. The cam becomes the critical element 
in the design; however, the accuracy requirement on the cam is not 
positional accuracy but that the incremental changes in the slope of 
the cam be no greater than 20$ of the average slope. This is true 
because the servo motor will finally drive the error detector package 
to the point where the stellar image impinges on the center of the 
knife edge. 

In Figure 64 a synchro is shown geared to the cam and servo motor. 

The purpose of this transducer is to indicate the satellite attitude error; 
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the attitude error in the plane 'which includes the synchro is ""he differ¬ 
ence "between the synchro signals from the two ends of the error detector 
package. 

Previously in this report, it was mentioned that the maximum displace¬ 
ment which the error detector would experience (normal to the optical axis) 
was plus or minus 0.10 inches, both in the horizontal and vertical direc¬ 
tions. This represents a capability of tracking stel?~r 'mages thru 
a range of 400 seconds of arc. Assuming tentatively that the cam will 
rotate 320° in order to provide the plus or minus 0.10 inch travel, if 
the synchro is geared directly to the cam, the scale factor at the synchro 
will be hJO seconds of arc of attitude error per synchro revolution. Con¬ 
sidering synchro accuracy alone, the capability then exists of measuring 
attitude errors to an accuracy of 0.1 seconds of arc. 

The motor-tach-gearhead unit and associated gearing, cam and synchro 
can be mounted as an integral sub-assembly. Two such sub-assemblies 
would nest into one area of the main aberrascope housing; an area which 
can be designed to be completely sealed from the ambient environment. 

The push-rods will of necessity be required to emerge from the sealed 
area; however, these rods can be sealed with mecal bellows since 
the range of linear motion is so small. Thus, there will be two sealed 
areas in the main aberrascope structure; these areas will be located at 
right angles to each other (with respect to the optical axis). 
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E. Photomultiplier Configurate t. 


27 

In the previous quarterly report, consideration was given to the 
arrangement of the photomultiplier tubes to provide a maximum of 
reliability. Two different types -,f mounting systems were considered. 

One was the conventional type :£ system which we called Type I where 
each stellar image was uniquely associated with only one photomultiplier 
tube. The incoming energy upon striking a knife edge prism would be 
directed to the photcmuitiplier tube after passing the u a chopper from 
either of two paths fr m the two faces of the knife edge. The two 
telescopic systems would be completely independent of each other as 
far as the detection system is concerned. Tf the probability of 
failure of one photomultiplier tube is designated as p , then the 
probability of failure of the entire aberrascope system (which involves 
two aberrascopes) expressed as is the following: 

* f - [p f (2 - P f )] 2 + 2 p f (2 - p f ) ( 1 - P f f 

In an alternate system which we called the l^ype II system, the stellar 

images from both stars were combined i” a unique chopping system such that 

each of the two photomultipriers receivea energy from both stars. While 

the Type I system blocked Out the stellar energy 50$ of the time, this 

(Type II) system utilizes all the incoming energy 100$ of the time. Again, 

if is the probability of failure of one photomultiplier tube, the 

probability of failure of the aberrascope system if is given as follows: 
h 2 , 2. 

7f t = P f + 2 p f ! 1 - P f 

The probability of failure functions for the two types of systems are 
shown in Figure 59 • The advantage cf the Type II system is shown quite 
forcibly in Figure 60. If the probability of failure of the individual 
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photomultiplier cubes is net greater than 0.1, then it can be seen that 
the Type II system has an advantage of a factor of 17 over the Type I 
system. It is probably doubtful whether a photomultiplier tube would be 
considered for use if its probability of failure was greater than 0.1 
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Probability of Failure - One Photomultiplier 
Figure 65 

RELIABILITY COMPARISON OF 
PHOTOMULTIPLIER ARRANGEMENTS 
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Figure 66 

RATIO OF PROBABILITIES OF FAILURE 
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F» Position Measuring Device 

It has teen mentioned already that the final output of the 
aherrascope will he the data obtained from measuring the displacement 
of the center of the error detector package? Thi* soasureaent will he 
accomplished ty means of a pair of interferometers: one in each of the 
two significant planes (horizontal and vertical). Hi-* interferometer 
will simply produce Newtonian fringes in a pattern of concentric circles. 

We ore assuming here that the central element on the error detector has 
s spherical surface so that the attitude changes of the satellite which 
will cause the error detector package to pivot about its center will not 
cause changes in the number of fringes. 

With the vise of helium light, the distance between fringes which is 
one-half wave length will he 11.6 x 10"^ inches. Tnen the change from a 
dark area to a light area will he one-half of this distance or 5*8 * 10“^ 
inches. By using a photo sensor to detect either a dark or light area, 
an electrical on-off type signal will be produced for every displacement 
larger than 5*8 x 10“^. In terms of angle in the aherrascope, this minimal 
displacement represents 0.012 saccods of arc. Since this measurement is 
made at the center of the knife edge, the actual error sensitivity will he 
twice ae large, or 0.02k seconds of arc. This illustrated in the sketch 
below. 


h- 1 —i 
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Thus, the error sensitivity capability is more than adequate to meet 
the requirements of the refraction measuring aberrascope. 

It can be recognized that a single photo sensor will not be 
able to distinguish between the increase or decrease of the distance 
between the fixed and moving elements of the interferometer. This 
problem can be circumvented by using two sensors: the spacing between 
positions will be such that while one is sensing a dark area the other 
is sensing a light area. The fact that the second sensor is observing 
either a light or dark area while the first sensor is changing from 
light to dark will establish the direction of the motion. 
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VI. CONFIGURATION OF STAS TRACKERS 


Since we are required to knew precisely which star is being 
occulted at a giver, time in order to derive useful information, it 
is necessary that we employ star trackers for collecting the basic 
data. By measuring both the time of ingress and the time of egress 
of each star it is possible to derive two transit* time measurements 
from each star tracker during each rotation cf the so- ll-te. Thus 
a minimum system requires that we use at least 3 star trackers. 

However, for certain orbits, 2 of the 3 stars can be occulted 
by the earth and s one star is inadequate to maintain the inertial 
orientation of the star trackers, a problem of reacquiring the stars 
is presented. If we add a fourth star tracker and arrange them as 
shown in Figure 67 at least two stars will be in the field of view 
at all times. Thus we will be able to maintain continuous altitude 
control and, depending on the orbit, we will have either h or 8 star 
transit time measurements per satellite, rotation. 

With the arrangement shown in Figure 67 the stars are picked in 
diametrically opposite pairs which are approximately orthogonal to 
one another. The star selection is made in this manner because of 
the possibility of mounting the two telescopes back-to-back and 
canceling their alignment errors by rotating them about their common 
optical axis; also in this way the refraction of the starlight which 
is passing through the earth's atmosphere can be more accurately 
measured if io is compared with its diametrically opposite companion. 


#Tne word "transit" in this caseis taken to mean "the passage of a smaller 
body across the disk of a larger". Thus an ingress and an egress of a 
given star represents two transits across a planetary limb. 


• I 


202 




2 STAR TRACKERS FOR 

CONTINUOUS ATTITUDE CONTROL 



Assume that the orbit is such that one palx- of stars is never 
occulted by the earth. In this case there will be only 4 star transit 
time measurements per rotation and we do not have enc ugh information 
to completely define the elements of the orbit. Since we can deter¬ 
mine the semi-major axis from one complete rotation, the remaining 
information enables us to calculate all the other orbital elements 
except the orbital inclination. Since the inclinatJ '0 vf the orbit 
relative to the equatorial plane is constant, a previous determina¬ 
tion of this parameter could he used. Another method of deriving 
the orbital elements in the special case where only one pair of stars 
is occulted is to measure the transit time of ingress or egress of 
the sun. This would involve the addition of some extremely simple 
instrumentation to the basic system and along with the star transit 
times would provide a basis for calculating all orbital parameters 
directly. 

In order to better understand the problem of star selection for 
various orbital planes, consider Figure 68. In this figure we have 
shown the transit boundaries for one, two and three aherrascopes, 

In drawing these transit boundaries a 300 n. mile circular orbit was 
assumed. In this case both members of each uiametric&ixy opposite 
pair can be tracked over a central arc o.f about 44 ° as shown. For 
the case in which one star pair is used, orbits which lie within 22* of 
the X-Y plane will involve no occultations. All other orbits will 
involve at least 4 star transits. For the case in which two star 
pairs are used, the least number of occultations is seen to he 4 . 

This occurs for orbital inclination lying within 22 c of either the 
X-Z o: Y-Z planes. For all other orbital inclinations 8 star 
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OF DIAMETR' CALL/ OPPOSED 30 O K.Ki. 
IrtCULAR 0R31T ALTITUDE STAR PAIRS 











transits will be observed. In the last case wixh 3 orthogonal star pairs 
either 8 or 12 star transits will be observed. The shaded areas of 
Figure 68 denote positions where no stars will be occulted. 

The transit boundaries shown in Figure 68 are dependent on the 
positions of the stars only and earth's pole or the plane of the orbit 
may be in any position relative to these patterns. For most orbits, it 
is only necessary to track 4 stars and for this cu*c ci : 4 or 8 star 
transit time measurements will be obtained per rotation. If the plane 
of the orbit should happen to intersect with the plane defined by the 
two pairs of diametrically opposite pairs, the error in the computation 
of the orbital plane, s\- and i will be very large as shown in Section IV-c 

The two major problems in tracking two diametrically opposite 
star pairs are that (1) in certain cases orbits can be obtained in 
which only 4 star transit time measurements are available, and (2) 
large errors will exist in the calculation of the plane of the orbit 
if it coincides with the plane of the stars. 
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In order ~o ciroumYen 4 ' these two problems one might consider a 
geometrical arrangement of 4 star trackers in which the angles between 
all of the stars are the same. 

',?ith the planar configuration shewn in Figure 69 as a starting condi¬ 
tion, assume o - 180% 0 = 90° • We now wish to increase")] until 0 = Q. 
Therefore, from the geometry we may write 

0 = 180 - 27 1 


Therefore 


cos (I) 


cosT) 

VT" 


-1 , CCS?] 

cos (-*- 


\nr 


180 - 27] 


tan 7? = —- 

~T\ = 35° 17’ 


Thus we now have an arrangement of star trackers in which the 
angles between all of the stars are the same and equal to 180° - 2(35° 17') 

= 109° 26'. This is shown in Figure JO . We now no longer have all 4 
stars lying in the same plane, but we have 6 different planes. Because 
of this there is no degeneracy in the calculation of any of the orbital 
elements. In addition, if the altitude of the orbit is not greater than 
800 n. miles, C star transits will be obtained for all orbital inclinations. 
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It is aot possible to specify a preferred, geometrical arrangement 
of the 4 star trackers which is the best one for a completely self- 
contained navigation. In general, one must consider: 

1) the semi-major axis, 

2) the eccentricity, 

3) the orbital inclination, 

4) the argument of the perigee, 

5) the longitude of the nodes, and 

6) the shape of the planet as it effects the rate of 
change of 4) and 5)* 

A configuration of star trackers should be selected which provides: 

1) a minimum of five stellar transits per rotation, 

2) the sinimum/maxiaum values of the error parametric 
sensitivity coefficients, and 

3 ) the simplest or most precise instruments for the 
measurement of star transit times. 
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